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HAUPTAUFSÄTZE 


Die elliptische Tragfläche auf potentialtheoretischer 
Grundlage’. 
Von Klaus Krienes in Göttingen. 


D‘ vorliegende Arbeit ist als Beitrag zur Theorie der tragenden Fläche gedacht. Es wird 
auf potentialtheoretischer Grundlage die Aerodynamik der elliptischen Tragfläche in ge- 
rader und schräger Anströmung untersucht. Die Grundlage der Rechnung bildet die lineari- 
sierte Theorie des Beschleunigungspotentials [1; 2]?), in der also alle Größen, die klein von 
höherer Ordnung sind, vernachlässigt werden. Das führt zu folgenden Vereinfachungen: 

1. Die z-Koordinate jedes Flügelpunktes wird vernachlässigt, d. h. die Unstetigkeit der 
Potentialfunktionen, die dem Drucksprung an der Fläche entspricht, liegt auf der 
Grundellipse (vgl. Bild 1). z 

2. Die Stromlinien, längs deren sd 
die konvektive Integration der 
Beschleunigung vorgenommen 


wird, sind Geraden parallel zur ; 
Anströmungsrichtung. 
Im Falle der elliptischen Beran- 


dung sind Lösungen der Laplaceschen 
Differentialgleichung Ay=0 als Pro- 


dukte von Lam&schen Funktionen be- 
kannt, deren Eigenschaften im folgen- 
den noch näher behandelt werden. Bild ı. 


Das Beschleunigungspotential und die Lameschen Funktionen. Die elliptische Tragfläche 
befinde sich in einer stationären Parallelströmung mit der Geschwindigkeit V. Die Flüssigkeit 
sei homogen, inkompressibel, reibungslos und der Schwere nicht unterworfen sowie außerhalb 
der tragenden Fläche und des abgehenden Wirbelbandes drehungsfrei. Es gilt dann 


DV 


1 


1) Dissertation Universität Göttingen 1939 (D 7). Berichterstatter: Prof. Dr. L. Prandtl. Der zweite Teil dieses 
Aufsatzes: „Der schiebende elliptische Tragflügel“ (S. 81ff.) wurde von der Lilienthal-Gesellschaft als Bewerbung zur 
Aufgabe 1 (Flugzeugbau) des Preisausschreibens 1938/39 angenommen und mit einem Preis ausgezeichnet. Eine aus- 
führlichere Fassung der hier veröffentlichten Arbeit befindet sich bei der Bibliothek der AVA Göttingen. 

2) Siehe das Schrifttumsverzeichnis am Schluß der Arbeit. 
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Der Druck p sei durch die Potentialfunktion y in folgender Weise gegeben: 


wobei p„ den statischen Druck im Unendlichen bedeutet. Die Funktion y genügt der La- 
placeschen Differentialgleichung 


2 


und man kann sie sich entstanden denken aus einer Quellsenkenbelegung der Tragfläche von 
der Stärke 0 (xp, yr,27r) an jedem Punkte der Fläche. 


/1 
Darin bedeutet » die Normalenrichtung der Fläche in (x, yr, 2r) und 


R=y(z — + (y— + 


die Entfernung des Aufpunktes (x, y,2) von dem Flächenpunkt (zr, Yr; 2r). Nach den anfangs 
gemachten Feststellungen kann das Integral ersetzt werden durch ein solches, das über die 
Fläche der Grundellipse genommen wird. 


Die rechnerische Behandlung des vorliegenden Falles wird ermöglicht durch die Ein- 
führung von ellipsoidisch- -hyperboloidischen Koordinaten [3]. Die Halbachsen der Grundellipse 
seien ce in der y-Richtung und cy1—x? in der «-Richtung (2xc ist dann der Abstand der 
Brennpunkte auf der y-Achse). Das neue Koordinatensystem wird nun eingeführt durch: 


2° 
x? 2? 


Die Auflösung der Gleichungen führt zu: 


c 
yzcouZ, 
(6). 
Ve 
Die Flächen o bzw. « bzw. »=const sind konfokale Flächen 2. Grades. o=1 liefert 


die elliptische Grundfläche, «=1 die Ebene z=0 außerhalb der elliptischen Scheibe. Das 
Flächenelement der Grundellipse ist 


dedy=e? 
Außerdem gilt für o=1 


Führt man für o, u, » mittels der Weierstraßschen ®-Funktion neue Koordinaten u, 
v, w in folgender Weise ein [4] 


indem man die in der Gleichung 


auftretenden Größen 
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d.h. 
setzt, so wird 
dp 


und die Laplacesche Gleichung lautet 


[P (e) —# (w)] (m) +[P (u) - 
Setzt man die Lösung in folgender Form an: 
-. . - 2.2. ....(i4, 
so erhält man für jede der drei Funktionen E die Lam&sche Differentialgleichung 
d’ E(u 
mit den Separationskonstanten A und B. Wir setzen jetzt 
B=n(n+1) und 


und führen mittels Gl. (9) wieder «u, v, o ein. Dann lautet die Lame&sche Differential- 


+ [1 vun —n(n+1) E,” 
Für B=n (n+1) gibt es genau 2» +1 Werte v,,, für die E„” („) die Form hat [6] 


. (15a). 


0 

E,” (u)= yi Vu — us (a, ur (16). 
0,” (u)=a,u" ist ein in gerades Polynom vom Grade n— 2,— Die 
Werte v„ sind voneinander verschiedene Lösungen einer algebraischen Gleichung, die sich 
aus der Bedingung ergibt, daß 0„”(u) ein Polynom ist. Mit Rücksicht auf die späteren 
Rechnungen sei E„”(»v) folgendermaßen definiert: 


wenn man 


v 


= (— cos p)* (sin . (16a), 


—=sing; _ ..... (MM 


setzt. Die mit Gl. (13) gewonnene Lösung 


läßt sich auch noch in anderer Form darstellen [7]. Bezeichnet man nämlich die Nullstellen 
des Polynoms 0,”(4) mit o,’, so sieht man leicht aus den Gln. (5) und (6), daß 


Aus der Klammer { } ist jeweils einer der darin enthaltenen Faktoren zu wählen, und 
das Produkt A ist über alle zu 0,”(zu) gehörigen Nullstellen zu bilden. Man kann die Null- 


stellen o,? sie direkt bestimmen, indem man auf die rechte Seite der Gl. (19) den Laplace- 


2 2 
schen Operator + anwendet und das Resultat gleich Null setzt. Das liefert 
folgendes für die 0,°: 
3 
—% Os 
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Die Potentialfunktionen in der Form (18) bzw. (19) sind sogenannte ‚innere‘ Lösungen 
der Potentialgleichung, die für unsere Zwecke nicht brauchbar sind, da wir Potentialfunktionen 
benötigen, die im Unendlichen durchweg verschwinden. Diese erhält man, wenn man in der 
Veränderlichen o nicht die Lamesche Funktion erster Art sondern die zweiter Art nimmt. 


Das ist diejenige Lösung der Lam&schen Gl. (15a), die für 0>» wie Sr >0 geht. Sie 
hat die Form 
do 
m m 
Das auftretende Integral läßt sich stets auf elliptische Integrale erster und zweiter 


Gattung allein zurückführen. Die äußere Lösung der Potentialgleichung sieht nun folgender- 
maßen aus: 


(21). 


(0, 1,7) = En” (u) Eu” (9) Fu" (0)= Eu” (1) Eu” (9) "oz, ©. 


Aus (19) wird dann 


st|, d 


zzx 0? 


Aus der Darstellung der Potentialfunktion durch (4a) als Quellsenkenbelegung auf der 
elliptischen Scheibe geht hervor, daß die Potentialfunktionen in der Ebene 2=0 im Außen- 
gebiet der Scheibe Null sein müssen. Nach (23) müssen also solche y„” genommen werden, 
die den Faktor 2 enthalten, d. h., von den Lameschen Funktionen kommen nur solche in 
Frage, die von der Form sind (vgl. (6) u. (16)): 


Die Lameschen Funktionen besitzen ähnliche Ollssnnlliäisigmnheftin wie die 
Kugelfunktionen [8]. 


+% 
d 0, t = 
+8 0,wennn=#+s 
od t 
und m=t 


- den Zusammenhang mit Gl. (4a) zu gewinnen, benutze ich folgende Darstellung 


2n+l 


R (u) Eu” (v) (0) Eu” (up) En” En” (op) . . (27), 


dann ist 


1 m m m m M„” (1) 


wobei dann nur noch über E„”(u„) von der Form (24) zu summieren ist. Nehme ich jetzt 
an, daß die Quellsenkenverteilung auf der elliptischen Scheibe am Rande mit der Wurzel 
aus dem Randabstand, d. h. mit YI— ur? gegen Null geht, so kann ich o (xp, yp) nach Pro- 
dukten von Lam 6schen-Funktionen der Art (24) entwickeln: 


Bilde ich jetzt das Integral (4a), so ergibt sich auf Grund der Orthogonalität der Lame- 
schen Produkte (vgl. (26)): 


n=l m 
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Gl. (29) zeigt, daß die Potentialfunktion (4a) eine bestimmte Linearkombination der 
Funktionen y„” ist, die zu untersuchen im folgenden die Aufgabe sein wird. 


Die Darstellung der Potentialfunktion ,„"” als bestimmtes Integral im Komplexen. Setzt man 


1 


(e, durch Gl. (11) definiert) oder mit Berücksichtigung von (6) und (9) 


so kann man zeigen, daß 


wenn man als Integrationsweg von x, y t-Ebene 
nach x, eine Schleife um die X=+1 um, m 
entsprechenden Punkte der komplexen t- 2 
Ebene nimmt. Q,„(X) genügt der Le- 
gendreschen Differentialgleichung 
+n(n +1) Q,(X)=0 
und E,„”(t) genügt der Lameschen 
Differentialgleichung (15). Es gilt dann ART, 
nämlich Bild 2. 
0° 0? Qu (X 
für je zwei beliebige der Variablen t,«,v,w, d. h. 
und 
1 X dE,„”(t 
Dies gilt auch, wenn ich « durch » oder w ersetze. Der Integrationsweg ist nun so 


ist, d. h. y„”" genügt in jeder der drei Variablen «, v, w der nn e schen Differentialgleichung (15) 
und ist demnach eine dritte Darstellung der durch (22) definierten Potentialfunktion. Die 
Punkte X=+1 in der t-Ebene sind gegeben durch 


t=v+n+u (X=—]); t=v+w—u (X=+H1l) .. . (8). 
Auf der elliptischen Scheibe (o=1) ist # (u)=e,, d.h. 
1 


do 


Als Beispiel sei nachfolgend die Potentialfunktion y,' angegeben. Die einzige mit dem 
Faktor Yo’ —1 versehene Lamesche Funktion erster Art und ersten Grades ist 


E'(0)= ve—i . 


3 
| 
| xum | 
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Nach (22) bzw. (23) lautet die Potentialfunktion dann 


d 


ı 1— x’ do 


Auftrieb und Auftriebsmomente. Der Auftrieb ist gegeben durch 


Dabei ist 
Yu |e=1= En” (u) Eu” (1) 


Auf Grund der Orthogonalität der l,ameschen Funktionen (vgl. 26) liefert nur y,' einen 
Beitrag zum Gesamtauftrieb 


y;' liefert das Kippmoment um die y-Achse 

y,* gibt das Rollmoment um die «-Achse 


(Der negative Index bezieht sich auf die in y ungeraden Funktionen.) 


Die gerade angeströmte elliptische Tragfläche. Die elliptische Tragfläche sei jetzt in 
Richtung der positiven x-Achse angeströmt mit der Geschwindigkeit V. Gl. (1) liefert nun 
die Möglichkeit, die von dem Druckpotential y im Raume und insbesondere an der tragenden 
Fläche induzierten Geschwindigkeiten zu berechnen. Die z-Komponente der Gl. (1) lautet im 
stationären Fall 


Größen, die klein von höherer Ordnung sind, werden vernachlässigt, d. h. 


Die z-Komponente w des Geschwindigkeitsvektors w soll im folgenden kurz mit Abwind 
bezeichnet werden. An der elliptischen Fläche erhalte ich den Abwind, wenn ich Gl. (46) 
für 2=0 und y<c über x integriere: 


Die Berechnung des Integrals wird auf elementarem Wege sehr mühsam und unüber- 
sichtlich bei größerem n. Wesentlich einfacher gelangt man zum Ziel, wenn man die Dar- 


m 
stellung (32) der Potentialfunktion y,„” benutzt. Bildet man ern durch Differentiation unter 
dem Integralzeichen, so wird aus (47) 


x 


2. . Math. h. 
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Nun ist 
 dX 92 HdX 

02 ' 
Setze ich dies in (48) ein und führe die Integration über x aus, sc erhalte ich, da 
lim Q,(X) =0 für ist: 


(49). 


2ri ( n ( ) 
Für Q@,(X) führe ich jetzt die Neumannsche Darstellung ein: 
+1 
[RAY 


Analog X(f) setze ich nun 


so daß ich für (50) erhalte: 


Der Integrand hat Pole an den Stellen t=s, da dort X(t)— Y(s)=0 ist und an der 
Stelle 2=w,+iw,, wo Y® (f)—e, eine einfache Nullstelle hat. Das Verhalten der Nenner in 
Nähe der Nullstelle wird bestimmt aus der Taylorentwicklung. Das liefert 


= 
unter Berücksichtigung von Gl. (10), sowie 
X Y 


nach Gl. (52). Führt man jetzt die Integration von s=v+mw-+tu bis s=ev+w— u ent- 
sprechend Y=—1 bis Y=-+1 (vgl. 38) durch, so wird 


Pu(Y ay 
Bei 
nm 


ist zu bemerken, daß der Integrand für o=1, d.h. u=—o, die Periode 2, hat, so daß 
ich den Integrationsweg so verschieben kann, daß s von iw,—o, bis io,+o, läuft (vgl. 


Bild 2), was bis +5 entspricht, wenn ich Ye ()—e,=—xsine, d.h. ds= 
(4e=yi —»’sin? e) setze. Weiterhin sei 
Es ergibt sich dann wegen 
(o,+io)=e) - - &9 


und nach (33 
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für die Abwindfunktion auf der elliptischen Scheibe 


Kir. 


Man erhält als Koeffizienten im Polynom von &£ und n vollständige elliptische Integrale 
1. und 2. Gattung. Die Berechnung der Abwindfunktion im Fall » = 1 sieht folgendermaßen aus: 


PNM=Y, 1 —1. 
Aus (61) ergibt sich 


en. =yi-:; ); E(3)=\4rae. 
Die tragende Fläche. Die Gestalt der Fläche sei durch 


gegeben. Die Neigung der Fläche in x-Richtung muß mit der Richtung der Strömung an 
demselben Punkt übereinstimmen, d. h. 


d2(2,9y)_ w(x,y) 


Daraus ergibt sich für 2=2z(«, y) 


x 


wobei die untere Integrationsgrenze beliebig ist; ich setze sie gleich Null und füge zum Wert 
des Integrals eine beliebige Funktion in y hinzu : 


x 


Die tragende Linie, der induzierte Widerstand und die Saugkraft. Der Kürze halber sei 
hier nur verwiesen auf die entsprechenden Abschnitte in der Arbeit von W. Kinner [2], in 
der diese Fragen ausführlich behandelt sind. Die Ergebnisse lassen sich ohne Schwierigkeit 
auf die hier behandelte elliptische Scheibe übertragen. Die tragende Linie, durch die man 
sich die tragende Fläche ersetzt denkt, erhält man, indem man die Auftriebselemente durch 
Integration parallel zur x-Achse zusammenfaßt: 


+XR 


Es ergibt sich dann z. B.: 
)=4r5 

| 
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Die Potentialfunktionen zweiter Art. Die im vorangehenden behandelten Potential- 
funktionen lassen sich in beliebiger Weise linear kombinieren. Man erhält dann für Auftrieb, 
Auftriebsmomente und den Abwind w die entsprechenden Linearkombinationen. Die bisher 
behandelten Potentialfunktionen liefern die aerodynamischen Größen einer entsprechend ge- 
wölbten Tragfläche bei stoßfreiem Eintritt der Strömung, d. h. bei einem bestimmten Anstell- 
winkel, bei dem keine Umströmung der Vorderkante eintritt. Den beliebigen Anstellwinkel 
erhält man, wenn man eine ebene elliptische Scheibe mit ihrer Strömung überlagert, bei der 
bekanntlich die Vorderkante Saugkante ist, d. h. die Auftriebsdichte geht dort gegen Un- 
endlich. Die Potentialfunktionen erster Art gehen jedoch alle am Scheibenrand mit der Wurzel 
aus dem Randabstand gegen Null. Es ist jetzt also die Aufgabe, Potentialfunktionen zu suchen, 
die die o. a. Eigenschaften besitzen. Diese erhält man, wenn man auf die Potentialfunktionen 
erster Art bei konstantem &, y, 2 und x? folgenden Grenzübergang anwendet (vgl. hierzu 
W. Kinner 


Diese Potentialfunktionen ®,„” haben die Eigenschaft, am ganzen Rand der Scheibe Un- 
endlich zu werden. Wegen der Bedingung des glatten Abflusses an der Hinterkante der Fläche 
ist es später notwendig, die Funktionen in entsprechender Weise linear zu kombinieren. 


Die Abwindfunktionen zweiter Art. Entsprechend Gl. (65) erhalte ich die Abwind- 


funktionen durch Anwendung des Öperators c" auf Gl. (61). Dabei ist Be- 
ziehung zu beachten 


dp dPu_, 
Dasselbe gilt auch für P„(Y) und Q,„(n). Aus (61) wird dann 
Po Ny 


„9 d 


day 


2 


Bei der Behandlung der zweiten Grundaufgabe, d. h. bei der Berechnung einer vor- 
gegebenen Tragfläche, braucht man Potentialfunktionen, deren Abwindfunktionen an der Scheibe 


von x unabhängig sind. Nach (46) bedeutet dies wegen a —=0, daß 


dx 
sein muß. Ich setze jetzt mit noch zu bestimmenden Koeffizienten 5,” 

m m 
Dementsprechend ist 
m 


w„(x,9) wird als Abwindfunktion der Potentialfunktion zweiter Art ®, bezeichnet. Nach 
(68a) und (66) ist 
+7 


cose 


Die Koeffizienten 5b,” sollen nun so bestimmt werden, daß w,„ eine Funktion von 7 
allein wird. Der erste Term in (69) hängt bereits nur von n ab; es muß also der zweite 


d Y 
Term, der im allgemeinen auch von £ abhängt, eine Funktion nur von n sein. rn) 


ist ein Polynom in Y mit Gliedern der Form 
ist also eine Summe von Gliedern der Form 
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Für die folgende Überlegung nehme ich an, M,„” (e) sei in e gerade; dann verschwinden 
alle Glieder mit ungeraden Potenzen von sin e bei der Integration von -5 bis +7. Ist a 


aber gerade, so ist 
(sin e)*=1-+ Summe von cos-Gliedern. 


Setze ich dies in (70) ein, so besteht also Y"-2» aus folgenden Summanden: 


mit der Bedingung, daß 
und verlange, daß 
+3 


so verschwinden alle Glieder, die Potenzen von cose, d. h. Potenzen von & enthalten. Den 
einzigen nicht verschwindenden Summanden von Y*-2» erhalte ich, wenn n —2q=0 ist, d.h. 
wenn erstens n gerade und nach (71) n —2p — a=0 ist; dieser lautet nach (70a) dann n%—?» 
und hängt nicht mehr von e ab, so daß ich ihn vor das Integral ziehen kann. Durch Zu- 
sammenfassung der Summanden vor dem Integral, das jetzt für jedes p denselben Wert hat, 


erhalte ich gemäß u AL in (69) für das Integral 


(Für n verschwindet das 

M„"(e) enthält den Faktor sine und ist demnach in e ungerade. Den vorangegangenen 

entsprechende Überlegungen führen dann zu der Bedingung 
+7 


| sine > =o, 
€ 


q=1,2...0-1; bw. 


Das Integral in (69) hat dann nur für ungerades » einen von Null verschiedenen Wert. 
Zusammenfassend erhält man unter Berücksichtigung von Ey,(x) = Ey," ı (x) = (vgl. Gl. (16)): 


1 1 
1 1 dQ 


Hierbei ist 
n 
+5 


1 
kır yı=— 2 +1 +1 | 


n 
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b„" und 5b„ ”" genügen den Bedingungen (73) bzw. (74). Das sind «—1 homogene 
Gleichungen für die o Unbekannten b„”, bzw. r—1 Gleichungen für die r Unbekannten b,", 
die daraus bis auf einen gemeinsamen konstanten Faktor bestimmt werden können. Letzteren 
will ich so wählen, daß 


wird. Dann ist 
1 #P,;- 
OF 1 (76). 


Auftrieb, Auftriebsmomente und die tragenden Linien der Potentialfunktionen zweiter Art. 


Diese Größen erhält man durch Anwendung des Operators Rn 2 c" auf die entsprechenden 


Größen der Potentialfunktionen erster Art (42, 43, 44, 64). Bei der vorzunehmenden Diffe- 
rentiation ist zu beachten, daß in dem Flächeninhalt der Ellipse Fy] der Faktor c? enthalten 


ist. Dann liefert unter Berücksichtigung von (68, 75, 75a) 


®, den Auftrieb A=85 
1 8 —P 
®_, das Rollmoment V? (79), 


E 3) 


Für die tragenden Linien ergibt sich: 


I 


a,()=425 V’eyi- Pin), 
2 =6, 

a,y)=0 


Die zweite Grundaufgabe der Tragflügeltheorie. Es handelt sich jetzt darum, die aero- 
dynamischen Größen eines beliebig vorgegebenen elliptischen Tragflügels zu berechnen ver- 
mittels der im vorangegangenen behandelten Potentialfunktionen erster und zweiter Art. Die 
zu erfüllende Randbedingung ist gegeben durch (62): 


%2(2,9) _ 
Differentiation dieser Gleichung führt zu (46) 
Setzt man y als Linearkombination von Potentialfunktionen erster und zweiter Art an 
v-2 
“= m n n 


so wird auf der Scheibe (o=1) 


02 


=0.. . .tla). 
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Nun ist aber 


En" (u) En") (0) 


Entwickelt man - in folgender Weise nach Lameschen Produkten 


so muß wegen Gl. (46), (81a) und (82) auf Grund der Orthogonalität der Lam&schen Produkte 


yi-—x? d(ye®— 1) Ay 
sein. Damit sind die Koeffizienten a,” der Potentialfunktionen erster Art bestimmt. Ent- 
wickelt man 2=2z(£,n) in eine Reihe nach & und n, so kommt die Bestimmung der a,” auf 
einen Koeffizientenvergleich in allen mit Potenzen von & behafteten Gliedern in Form eines 
Siebverfahrens hinaus (vgl. W. Kinner [2]. Wegen der Berechnung der a,” aus Gl. (46) 
ist nun 


eine Funktion von n allein. Diese Restbedingung wird nun mittels einer geeigneten Linear- 
kombination der Abwindfunktionen zweiter Art erfüllt, die ja nur von n abhängen. Die Be- 
dingung lautet: 


Diese Gleichung wird von „=0 bis n integriert und dann — mit P,.-, (n) bzw. Ps, (n) 
multipliziert — nochmals von „= — 1 bis +1 integriert. Das liefert unter Benutzung folgender 
Beziehungen: 


= wenn r=a 

0, wenn r=Fa 
+1 


1 


zwei unendliche Gleichungssysteme für die Koeffizienten C„ und D, 


\ 
y=1,2,... 


Die Abflußbedingung. Neben der Erfüllung der Gleichungssysteme (86) müssen die 
Koeffizienten C„ und D, noch der Abflußbedingung genügen, d. h. sie müssen so gewählt 
werden, daß die Auftriebsdichte an der Hinterkante verschwindet. Für die Potentialfunktionen 
erster Art ist diese Bedingung ohnehin erfüllt, da diese am ganzen Rand verschwinden. Das 
Verhalten der ®„” am Rand ist gegeben durch Anwendung des Differentiationsprozesses 
(65) auf (41): 

uzl c — 


I-»’ 

(Der Ausdruck ((y1— u?)) verschwindet mit der Wurzel aus dem Randabstand.) Ich setze 
m m 
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Dann wird für die Potentialfunktion zweiter Art 


Die Funktionen M, (9) besitzen nun ähnliche Orthogonalitätseigenschaften wie die tri- 
gonometrischen Funktionen, wie sich auf Grund der Bedingungen (73) und (74) nachweisen 
läßt, wenn man berücksichtigt, daß sich M,(9) in der Form 

M,(p)= (sin [a, (cos + a, (cos pP +...] 

schreiben läßt. Es ist 

# 

wenn n—y=+23,+4... .... (89. 

Die Abflußbedingung lautet nach (88): 


r=1 

3 
(+2 Ders M-er =0 (0b). 


Diese Gleichungen multipliziere ich mit (9) Ay bzw. M_2,(g) und integriere 


von D. h. 


Day Der I-er y=0 v-12,... . (92b). 
Dabei ist 
+ — 
2 


Die Systeme der unendlich vielen Gl. (86) sowie (92) liefern jetzt die Möglichkeit, die 
Koeffizienten der Potentialfunktionen zweiter Art C„ und D,„ zu bestimmen. Mittels dieser 
und der bereits berechneten Koeffizienten a,” der Funktionen erster Art ist man dann in 
der Lage, die aerodynamischen Größen der vorgegebenen Tragfläche zu berechnen. An einem 
Beispiel wird der Rechengang erläutert. ® 


Bemerküng. Die Rechnung des vorliegenden Teils ist durchgeführt unter der Annahme, 
daß die elliptische Scheibe parallel zur kleinen Hauptachse angeströmt wird. Die Länge der 
Hauptachse in y-Richtung war ce, in «-Richtung eyi—x?. Man kann die Rechnung ohne 
Schwierigkeit auf den Fall übertragen, daß die elliptische Tragfläche parallel zur großen 
Hauptachse angeströmt wird, indem man x? negativ annimmt, d. h. den Brennpunktabstand 
auf der y-Achse imaginär ansetzt. Für die elliptischen Integrale mit negativen x? gelten 
noch folgende Umrechnungsformeln: 


IR 


2 


v 


Die Integrale sind, bis auf den Faktor vor dem Integral, dieselben wie beim reziproken 
Achsenverhältnis. 
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Die ebene elliptische Scheibe in gerader Anströmung. Der Anstellwinkel sei «,, d.h. die 
elliptische Fläche ist gegeben durch 


— rg, = 

Nach (62) ist dann 
de. 


auf der Scheibe. Nach (83) sind deshalb die Koeffizienten a,” der Potentialfunktionen erster 
Art sämtlich Null. Die Potentialfunktionen zweiter Art müssen jetzt gemäß den Bedingungen 
(73; 74; 75a) bestimmt werden, was hier der Kürze halber nur für n=3 und n=—3 gezeigt 
werden soll. 

Das Achsenverhältnis sei 


vi-#"=- =0.%6, d.h = 6.37 
Für die vollständigen elliptischen Integrale erhalte ich dann 
# 
Nach (16a) ist im Falln=3 und 8, =1; ,=.,—=0 
s=1;2, — 0°, (e)=x* sin’ e — 
1 1 3 
Nach (20) ist + =0, d.h 0, —=0.97008 . 


Bedingung (73) und (75a) lautet jetzt: 


ja 
a 


3 3 
b | (0.96 sin? — 0.97008) Ae de 0 


— — (0.96 — 0.19792) — — (0.96 — 0.97008) . 


Daraus ergibt sich 
b, = 0.309 . 


D.h. 
M, 1.627 sin? — 0.361. 
Für n=3 und 8=,=8,=1 erhält man 


Bedingung (75a) lautet: 


nt 
de 
cose 
d.h. 
b,'Y=2654. 
Demnach ist 
M_,(p) = — 2.654 cos psinp, 


In analoger Weise geht die Bestimmung der weiterhin benötigten Funktionen vor sich. 
Mittels der Funktionen M,(g) können nun die Integrale (93) berechnet werden. Diese sind 
entweder elliptisch oder lassen sich auf elementare zurückführen. Die Koeffizienten I/,,; der 
Gleichungssysteme (92) sind damit bekannt. Die Berechnung der rechten Seite von (86) ergibt 


0 für a>1. 
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Für y=1,2,... ist die rechte Seite von (86) stets Null, so daß in y unsymmetrische 
Potentialfunktionen nicht auftreten. In der durchgeführten Rechnung wurde die Reihe (81) 
der Potentialfunktionen mit n=4 abgebrochen, so daß die Gl. (%) und (86) für a=1;2 
genommen werden müssen. Das liefert 


2.101 C,+ 1.5217 C, +0527 0,=0 (@=1)} 
— 0.2410 C, + 0.5132 C, + 0.9248 0, —0 
1 2 1 2 
(86). 
15 (a=2) 
Die direkte Auflösung des Gleichungssystems führt zu 
C,=0.568 a,, C,= — 0.342 a,, C,=—0.0135 o,. 
Nach (77) ist jetzt der Auftrieb 
A=0568a,:8:5 = 455,5 V° Fan, 
d.h. 
da, > 4.55 . 
Für das Moment um die y-Achse ergibt sich nach (78) 
—= 0.7785 a,: 0.9524 ceyi— x? S V? Fan, M=—1%8a,cyl— x»? V? Fon: 


1—x 
Zur Berechnung des induzierten Widerstandes muß bemerkt werden, daß die Auftriebs- 
verteilung der tragenden Linie (vgl. S. 75) bei Berücksichtigung von nur endlich vielen Reihen- 
gliedern in (81) an den Flügelenden nicht verschwindet (vgl. Bild 3), so daß für die Be- 
rechnung eine Ersatz- Auftriebsverteilung 


Der Druckpunkt liegt bei vim= —0.435, d.h. in 28.3° max. Flügeltiefe. 


08 
genommen werden muß. In diesem Fall En 
eignet sich am besten hierfür die ellip- re nn 
tische, so daß sich für den induzierten 06 I 
Widerstand lediglich die bekannte Formel 
ergibt: errechnefe Auftriebsverteilung 
7 ptische Auftriebsverteilung — 
N 
Für das Achsenverhältnis | | ar 


#=05, 4-25 


gehen die Rechnungen in derselben Weise vor sich, indem man zunächst die Größen b„” aus 
den Bedingungen (73; 74; 75a) bestimmt. Damit hat man die Funktionen M, (9) für die 
Berechnung der Integrale I,,;, womit die Koeffizienten der Gleichungssysteme (9%) gegeben 
sind. Die Auflösung liefert 


C, =0.3741a,, C,= — 0.6347 a,, C,= —0.2347 a,, — 0.0138 a,. 
Der Auftrieb ist dann 
dca 


o 


Für das Kippmoment ergibt sich 
=— 1397 Fan. 
Der Druckpunkt liegt bei 


— 0.467 gleich 26,7°/ max. Flügeltiefe. 
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Für das Achsenverhältnis 
yi-«'=2 —3 A = 0.637 


ergibt sich 
C‚=0.124a,, 0,=—05245a,, 0, = — 0.066 a,, C,= —0.0lla,, 
d.h. 
0.9. 


Der Druckpunkt liegt bei 
0584 gleich Tiefe. 


Zum Vergleich seien noch die Werte für die ebene Kreisscheibe angegeben (vgl. 
W. Kinner [2]): 


»=0, A=1272, -18. 
Druckpunkt: = 0515, d.h. in 24,3°% Tiefe. 


Die in dieser Arbeit angewandte Berechnungsmethode läßt auch einen Grenzübergang 


zur tragenden Linie «”=y1-—x°’=0) zu. Man erhält dann bei Berücksichtigung von zwei 
Reihengliedern 


Druckpunkt in 28,8°) Tiefe. 
o 
Letzteres Ergebnis entspricht einer elliptischen Auftriebsverteilung über der Spannweite mit 


einer Druckpunktlage in + Tiefe in jedem Flügelschnitt (Schwerpunkt einer homogenen Halb- 


ellipse). Entwickelt man alle Größen, die in den Gleichungssystemen auftreten, nach x’, so 
erhält man folgende Näherungsformel für kleine x’: 


den _ 
da, 


Berechnet man ge nach der Theorie der tragenden Linie, nach der bekanntlich 
o 
1 dc, _2nA 


ist, so zeigt sich bei kleinem A eine erhebliche Abweichung gegenüber den nach der Theorie 
der tragenden Fläche berechneten Werten, 


da vgl. hierzu Bild 4. 
ze Andererseits ist die Übereinstimmung mit 
60 den von Weinig [10] errechneten Werten gut. 
2 \ In Bild 5 ist die Abhängigkeit der 
Druckpunktlage vom Seitenverhältnis dar- 
gestellt. 
N 
| 
A 
l A 
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Der schiebende elliptische Tragflügel. 
Unter denselben Voraussetzungen, die - 
bei gerader Anströmung der Tragfläche 
gemacht wurden, läßt sich auch die N 
schiebende elliptische Tragfläche be- 
handeln. Der sogenannte Schiebewinkel / 
sei 8, dessen Definition sich aus Bild 6 
ergibt. 

Die Stromlinien sind jetzt Geraden 
in der Ebene 2=0, die durch gr 


y=—(z— ap)tgß+y 7 
= — xtgß-+ const 


definiert sind. Aus Gl. (1) 


Wegen (%) ist aber 


Gl. (96) läßt sich u in der ey schreiben: 


Den Abwind erhalte ich daraus durch Integration längs einer Stromlinie, die wieder als 
Gerade parallel zur Anströmungsrichtung angenommen wird. 


12 cos 


Die Potentialfunktion y wird wieder als Linearkombination der Potentialfunktionen 
Yu", ®„,®_„ in der Form (81) angenommen. Für Auftrieb und Momente erhalte ich also 
dieselben Formeln. 

Für die Berechnung der Abwindfunktionen führe ich in der x y-Ebene ein neues schief- 
winkliges Koordinatensystem ein, das durch den zur Anströmungsrichtung parallelen Ellipsen- 
durchmesser — &;Achse — und den dazugehörigen konjugierten Durchmesser — ng.Achse — 
gegeben ist. Setze ich 


— x? si 
=yil—rtgß, d.h. sin (38), 
so sind die eben erwähnten Durchmesser gegeben durch 
yi-» . . . . . . . 
Auf der Scheibe, d.h. für o=1, ist nach Gl. (6) und (17) 
x 2 — x? 


In analoger Weise setze ich 
cos(P — 73 = usin (pP — pp), 
d.h 


csß, yi-—x»’sinß 


x? 


sinß, , cos ß 
Aß Aß n . . . . (100). 


| 

INeosß 

Bild 6. 
wird unter den vereinfachenden Annahmen 
©; 
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Ich benutze jetzt wieder die Darstellung (32) der Potentialfunktionen y„”. Bei den 
nachfolgenden Differentiationen bzw. Integrationen längs einer Stromlinie ist stets zu beachten, 
daß dabei y= — + const ist. 


=, 


2 
Ve,—e, 


(32). 


Es ist 
AQ„(X) 
02 
und 

ar (X) läßt sich also wieder durch — —— _ = ausdrücken; setze ich das Resultat in (97) ein und 
führe die Integration über x aus, so Be sich 


‚ts? 


Unter Benutzung von Gl. (51) und (52) wird wieder Q,(X) ersetzt durch 


Yıs)ds 


v+w— 


ds. 


Der Integrand hat Pole an der Stelle t=s, da dort X (t) —Y(s)=0 ist, und bei t=1ß, wo 


ps; Velt)—e, = — Ve,—e,c08 93. 


Der Nenner verhält sich an dieser Stelle wie 


ist, d.h. 


Das Residuum an der Stelle ist also 


os 


icospatitg ps sing; En 
unter Berücksichtigung der Tatsache, daß nach Gl. (100) 


X (t) = sin 


ist und nach (33) 


Für das Residuum an der Stelle t=s gilt wieder das bereits im Fall der geraden An- 
strömung Gesagte (vgl. S. 71). In diesem Fall erhält man dafür 


+5 
ide 
cose— Yl—»’tgßsine 


Ersetze ich noch Yi—x’tg £ durch tg 9; und führe mittels Gl. (100) die Koordinaten & und ns 
ein, so wird mit den Werten (102) und (103) aus Gl. (101) 
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Bd. 20 Nr. 2 April 1940 


1 
En” (— x 008.95) On (np) — 


(109), 


Bei den Gliedern von P,„, die Potenzen von £; enthalten, hebt sich der Nenner cos (e+ 95) 


im Integral fort, und man erhält vollständige elliptische Integrale erster und zweiter Gattung. 
Für die Berechnung der übrigen erweitere ich Zähler und Nenner mit cos(e — 95) 4? und 


erhalte dann mit 
cos (£ + 93) cos (e — 9) 4? B = cos’ B — 4? Bsin? e 


1 
AB E„” (— x c08 93) Qn — 
(104). 
1 
\ cos (+93) + ngsin (e + B dede 


Die auftretenden Integrale lassen sich zurückführen auf vollständige elliptische Integrale 
erster und zweiter Gattung und ein Integral der Form 


2 
sin? de 
cos? — 4? Bsin’e de’ 


= 


das sich als vollständiges elliptisches Integral dritter Gattung auf unvollständige erster und 
zweiter Gattung zurückführen läßt [11]. 


(os. 


a 


cos ß 
Für die Potentialfunktion y,' soll jetzt die Berechnung der Abwindfunktion als Beispiel 


durchgeführt werden. 


X, X+Hl 


x 008 9))= x? cos’ 

Dann ist nach (104a) 


cos (e — 95) de 


1 
v 
— 
+5 +, 


sin sin cos e 
cos? — 4? sin? e 

n 


= 
| 
. 
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Die Auswertung der N liefert 

1 

Die Abwindfunktionen der Potentialfunktionen zweiter Art. Nach (68) ist &,= Ib," ®,”, 
wobei nach (65) 

d 
Du" de [e” 9,2, e)] 

ist. In derselben Weise berechnet sich die Abwindfunktion w, aus w„” (&3,n3), 
M, (e) 


cos (+93) A ede (106). 


0, auf einer Stromlinie = const) ist die Abwind- 


funktion konstant. w, ist also eine Funktion von ns allein. Ich kann demnach in ——; P„_,(Y) 


dp 
dY 
&s—=0 setzen. Zu untersuchen. bleibt noch der Summand 


Setze ich wieder ı=5 bzw. ei, so wird 


[sin (e + = [sin (e + =?*- [sin (e + 
= [sin (e + ,)]" -2”-(1 + Summe von cos-Gliedern). 


Aus der Summe kann ich den Faktor cos (e+ 9;) einmal herausziehen, so daß er sich 
gegen den Nenner kürzt. Die neu entstehende Summe, multipliziert mit [sin (e + 9)" -?", 
liefert aber nur Glieder, die entweder in e ungerade oder der Bedingung (73) bzw. (74) 
genügen, und demnach alle verschwinden. Einen von Null verschiedenen Beitrag liefert nur 
[sin (e +93)]" —-?"-1, das aber vom Exponenten p nicht mehr abhängt. Ich kann also alle 


n3"- vor dem Integral gemäß (Y) zu zusammenfassen, da das Integral 
für jedes p dasselbe ist, d.h. 


n 1 m 
= — Aß b„ Eu" (— x cos + 


+, (107). 
n 


2 
Die Abwindfunktionen bei der schräg angeströmten Tragfläche haben demnach folgende 


Form: 


Die zweite Grundaufgabe im Fall des schiebenden Tragflügels. Man geht vor wie im Fall 
des gerade angeströmten Tragflügels, indem man 


dz(x, y) 


setzt, wobei 


y) = cos B . (108a) 


m n 
Auf der Scheibe ist 


Math. 
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und y=— xtgß + const ist. y sei wieder durch 
v= Ya," y" +3 +28 (109) 


gegeben. Gl. (96a) liefert dann 
d (e) 


womit die Koeffizienten der Potentialfunktionen erster Art bestimmt sind. Für die Be- 
stimmung der Koeffizienten der Funktionen zweiter Art nehme ich Bedingung (97) in folgender 


Form: 


Die rechte Seite ist eine Funktion von 7; allein wegen der Bestimmung der a,” aus (%a). 
Die Restbedingung lautet dann also: 


d n-ı n-ı 


Ich integriere diese Gleichung von 73=0 bis 7; und dann — mit P3a—ı (ns) bzw. Pa, (ns) 
multipliziert — nochmals von —1bis+1. Auf Grund der Örthogonalitätseigenschaften 
der Kugelfunktionen ergibt sich dann 


(1). 


(112). 


(113a). 
—10 
—10 


Die Abflußbedingung. Für die Formulierung der Abflußbedingung sei folgende Be- 
merkung eingeschoben. Die Abwindfunktionen (107) beim schiebenden Tragflügel haben die 


Eigenschaft, bei wie Q1n-ı(ns) bzw. unendlich zu werden. n3—=+1 sind 


aber die Randpunkte der elliptischen Scheibe, in denen Parallelen zur Anströmungsrichtung 
die Ellipse berühren (vgl. Bild 6). Man wird demnach anzunehmen haben, daß die soge- 
genannten ‚Wirbelzöpfe‘ in diesen Punkten parallel zur Anströmungsrichtung von der Scheibe 
abgehen. Ich verlange deshalb, daß zwischen den Punkten n3=+1 keine Umströmung der 
Hinterkante stattfindet. D. h. (vgl. S. 76): 


3 
te. 


Diese Gleichung multipliziere ich mit 
Ap bzw. M_2,(p) Ap 


und integriere von 3 +93 bis hm. Das liefert folgende Gleichungssysteme: 


86 Krienes, Die elliptische Tragfläche auf potentialtheoretischer Grundlage Az Ar 


r=ı 


(114). 
r= 
Dabei ist 
. (Ilka). 


Die Gl. (113) und (114) liefern jetzt die Möglichkeit, die Koeffizienten C, und D,„ zu 
bestimmen. Mit den bereits oben bestimmten a,” ist man dann in der Lage, die aero- 
dynamischen Größen einer schräg angeblasenen elliptischen Tragfläche zu berechnen. Zu den 
Gleichungssystemen ist noch zu bemerken, daß sie eine Kopplung der Systeme (86) und (92) 
für den Fall der geraden Anströmung bilden. Man erkennt dies leicht an der Ähnlichkeit der 
entsprechenden Gleichungskoeffizienten, womit man eine erste einfache Rechenkontrolle hat. 


Die schräg angeblasene ebene elliptische Tragfläche. Die Rechnungen sind durchgeführt 
worden für das Achsenverhältnis y Iw=+ (vgl. S. 78ff.) und die Schiebewinkel 


ß=15° und B=W0°. 

Benötigt werden hierfür wieder nur die Potentialfunktionen zweiter Art gemäß (108) und 
(110). Zwecks Vereinfachung der sehr umfangreichen Rechnung werden nur die Funktionen 
bis zum Grad n=3 genommen, d. h. 

D,, D_, D_, (®_, existiert nicht). 
Für die 5 unbekannten Koeffizienten dieser Potentialfunktionen brauche ich also 


5 Gleichungen. 2 nehme ich aus der Bedingung zw (n)=— eos ß tg a, auf der Scheibe, 


Gl. (113a) für a=1 und (113b) für y=1. Für die drei weiteren Gleichungen nehme ich 
die Abflußbedingung (Gl. (114) mit «=1,2 und y=1). In die Berechnung der schiebenden 
Fläche gehen dieselben M,(P) ein, die schon im Fall der geraden Anströmung bestimmt 


wurden. Die Integrale 10, erhalten jetzt nur die Grenzen 5+ 9; und top . Die Inte- 


grale I,,, bleiben dabei wegen der Periodizität des Integranden unverändert. Die für die 
Bestimmung der Koeffizienten k,„(ß), i„(#) usw. notwendigen Werte der unvollständigen 
elliptischen Integrale erster und zweiter Gattung wurden den Legendreschen Tafeln [12] 
entnommen. Für $=15° ergibt sich dann 


C,= 0.5204 a, 
(,= — 0.1719 a, D,= 0014 «, 
C, = — 0.3343 a, D, = — 0.0065 a,. 


A= 116,5 


M=—183a,ce 


2 
L= 00384 V* 
Für = 30% ist 
04070, 
C,= — 0562 a, D,= 0.0277 a, 
02650, D,= — 0.0124 a,. 


A= 3.260, V? 


5 Fa 
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Das neu hinzutretende Moment L ist nach den obigen Berechnungen positiv, was damit 
gleichbedeutend ist, daß die voraneilende Flügelhälfte größeren Auftrieb erhält. Nachfolgend 
sind noch die Koordinaten der Druckpunkte angegeben: 


2=0005  (p= 6°) 
0489; 70.097 (P=13° 


Zum Vergleich seien an dieser Stelle noch einmal die Werte angegeben, die man erhält 
bei #=0*°, wenn die Reihe der Potentialfunktionen mit na —=3 abgebrochen wird. In Klammern 
ist die Änderung in Prozenten gegenüber den mit 4 Reihengliedern berechneten Größen an- 
gegeben, woraus sich Schlüsse bezüglich der Konvergenz der Rechnung ziehen lassen. 


C= 056a, (-09°,) 


(-02°],) D,—=0 
(47.09) 
D. h 
A= 45005 
5 
Druckpunkt: = 0438 (+0.7°],). 


Eine Zusammenstellung der Ergebnisse zeigen die Bilder 7 und 8. 


d 
. keit vom Schiebewinkel 3 für das | 
Seitenverhältnis A = 6.37. 
[1-6] Bild $ (rechts). dee Abhängig- | 
da | | 
Seitenverhältnis 4 = 6.37. | | | 


Für große A läßt sich wieder eine Näherungsformel für das Rollmoment in Abhängig- 
keit vom Schiebewinkel $ und dem Achsenverhältnis x’=y1-—x? angeben. 


1 cos 128 * 


Die Formel zeigt, daß das Moment beim Grenzübergang zur tragenden Linie (x’ = 0) 
verschwindet. Macht man jedoch das Moment mit der halben Flügeltiefe an Stelle der 
halben Spannweite dimensionslos, so daß der Faktor x’ in (115) wegfällt, so wird der 
Momentenkoeffizient beim Grenzübergang logarithmisch unendlich. Im Grenzfall nimmt der 
Auftrieb mit cos’ ö ab, was sich dadurch erklären läßt, daß die a in 
«-Richtung V cos ß ist. 
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Die von S. Hoerner [13] veröffentlichten Messungen an schiebenden Flächen mit ver- 
schiedenen Umrißformen stimmen sehr gut mit den hier berechneten Werten überein. Es 
erscheint daher die Annahme berechtigt zu sein, daß das Rollmoment neben dem Schiebe- 
winkel im wesentlichen von der Flügelstreckung und nicht von der Tiefenverteilung im 
einzelnen abhängt. 


Zusammenfassung. In der vorliegenden Arbeit ist die elliptische Tragfläche in gerader 
und schräger Anströmung auf potentialtheoretischer Grundlage behandelt worden. Ent- 
sprechend der Theorie erster Näherung, die die Grundlage der Rechnung bildet, müssen die 
bekannten Forderungen über die Gestalt der Fläche und deren Anstellwinkel erfüllt sein. 
Als weitere Bedingung kommt noch hinzu, daß die Neigung der Fläche in Richtung der 
Stromlinien eine stetig differenzierbare Funktion der Flächenpunkte sein muß. Ist dies in 
einem vorgegebenen Beispiel nicht der Fall (z. B. bei Querruderausschlag oder V-Form der 
Fläche — Letzteres ist bei der schiebenden Fläche von Bedeutung —), so muß man sich be- 
helfen, indem man die Knicke durch entsprechende Abrundungen ersetzt. Im allgemeinen 
Fall wird für die Berechnung einer vorgegebenen elliptischen Fläche eine Reihe von unend- 
lich vielen Potentialfunktionen benötigt, deren Koeffizienten sich aus linearen unendlichen 
Gleichungssystemen bestimmen. Um zu einem zahlenmäßigen Ergebnis zu gelangen, bricht 
man diese Reihe gemäß der gewünschten Genauigkeit mit einem bestimmten Glied ab. Auf 
. die Integralgrößen Auftrieb und Auftriebsmomente hat dies einen praktisch zu vernach- 
lässigenden Einfluß. Über den induzierten Widerstand läßt sich keine unmittelbare Aussage 
gewinnen, da man hierfür sämtliche Koeffizienten der unendlich vielen Potentialfunktionen 
brauchte. Die Auftriebsverteilung geht nämlich sonst an den Flügelenden rechnerisch nicht 
gegen Null oder der Abwind wird dort unendlich, was dadurch bedingt ist, daß z. B. die 
Auftriebsverteilung der tragenden Linie hier nach Kugelfunktionen (vgl. (80)) entwickelt;wird, 
die an den Flügelenden nicht zu Null werden wie die sonst verwendeten trigonometrischen 
Funktionen. Bei der schiebenden Fläche, die man nicht ohne weiteres durch eine tragende Linie 
ersetzen kann, müßte man zunächst den sogenannten Flächenwiderstand (vgl. W. Kinner [2]) 


Wr=— de ay 


bestimmen und dann davon die Saugkraft an der Vorderkante subtrahieren. Dabei ist man 
jedoch wieder wegen der endlichen Anzahl der berechneten Koeffizienten auf Extrapolationen 
angewiesen. Auch die sich ergebende Druckverteilung ist bei wenigen Reihengliedern nur 
sehr bedingt gültig, insbesondere in Nähe der Flügelenden. 


Erwähnt sei noch, daß die in dieser Arbeit berechneten Potential- und Abwindfunktionen 
beim Grenzübergang x’>0 in die von Kinner bei der kreisförmigen Tragfläche verwendeten 
entsprechenden Funktionen übergehen. 


Die zum Teil umfangreichen Zahlenrechnungen wurden größtenteils mit der Rechen- 
maschine durchgeführt, da bei der Berechnung der elliptischen Integrale für höheres » Diffe- 
renzen von großen Zahlen auftraten, so daß die Rechenschiebergenauigkeit nicht mehr aus- 
reichte. 127 
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Elastizitätstheorie zweiter Ordnung mit Anwendungen. 
Von Maurice A. Biot in New York. 


Einleitung. In der klassischen Elastizitätstheorie geht man von der Voraussetzung aus, 
daß die Beanspruchungen, Formänderungen und Drehungen sehr klein sind. Sie werden 
durch Glieder erster Ordnung dargestellt, während Glieder zweiter Ordnung vernachlässigt 
werden. Im folgenden wird eine Theorie entwickelt, die alle Glieder zweiter Ordnung enthält, 
jedoch nicht die explizite Formulierung der Beziehung zwischen Beanspruchung und Form- 
änderung erfordert. 

Die Resultate wurden mit Hilfe der Variationsrechnung abgeleitet, jedoch ist das ganz 
unwesentlich. Der Erfolg dieser Methode ist darauf zurückzuführen, daß für die Komponenten 
der Formänderung neue Ausdrücke (Abs. 1) gewählt wurden. Dadurch ergibt diese Methode 
Glieder, die sich in der klassischen Theorie nicht vorfinden, die man jedoch mit Hilfe der 
Gleichgewichtsbedingungen eines kleinen Elementes auch erhalten kann. Die Komponenten 
der Formänderung stehen in linearer Beziehung zu den tatsächlichen Längenänderungen. 


Damit wird ein Ausdruck für die Dehnungsenergie aufgestellt (Abs. 2). Es ergibt sich 
in natürlicher Weise die Dualität der Darstellung des Beanspruchungszustandes. Zugleich 
werden Beziehungen abgeleitet zwischen den Komponenten der Beanspruchung, bezogen auf 
die Fläche vor der Formänderung, sowie denen, bezogen auf die Fläche nach der Form- 
änderung. 

Im Abschnitt 3 sind die Gleichgewichtsbedingungen abgeleitet. Diese enthalten alle 
Glieder zweiter Ordnung und beziehen sich auf beide Systeme der Beanspruchungen: Auch 
hier wurde die Variationsrechnung angewandt. 

Im letzten Abschnitt wird die lineare Theorie für kleine Formänderungen und kleine 
Änderungen der Beanspruchungen für einen Körper unter anfänglicher Beanspruchung be- 
handelt. Daselbst wird ein Ausdruck für die potentielle Energie mit allen Gliedern zweiter 
Ordnung hergeleitet, sowie die linearen Gleichungen für das Gleichgewicht. 


Diese Theorie ist vor allem anwendbar bei Behandlung der Fortleitung von Schwin- 
gungen in einem vorbeanspruchten Körper und bei Untersuchung des elastischen labilen 
Gleichgewichtes und der Knickung. Auch hier müssen wir streng unterscheiden zwischen 
kleinen Änderungen der Beanspruchung, bezogen auf die Flächen vor oder nach der Form- 
änderung. Es ergab sich, daß die Elastizitätsmoduli des Hookeschen Gesetzes für die zu- 
sätzlichen Beanspruchungen und Formänderungen eine symmetrische Matrix sind, wenn die 
Beanspruchungen auf die Flächen vor der Formänderung bezogen sind. Dieses ist nicht der 
Fall für Beanspruchungen, bezogen auf die Flächen nach der Formänderung. Für den vor- 
beanspruchten Körper ergeben sich Gleichgewichtsbedingungen, die, für Beanspruchungen 
bezogen auf die Flächen vor der Formänderung, schon von C. Biezeno und H. Hencky') 
auf ganz anderem Wege erhalten worden sind. 

Nahe verwandt mit dieser Aufgabe ist das Problem der elastischen Stabilität, das auch 
schon verschiedentlich untersucht worden ist. R. V. Southwell?) hat den Fall einer an- 
fänglichen gleichmäßigen Vorbeanspruchung behandelt mit den Koordinatenrichtungen längs 
der Hauptachsen. E. Trefftz*) hat Gleichungen aus der Dehnungsenergie abgeleitet, die 
den genauen Ausdrücken nicht ganz entsprechen. Erwähnt seien auch die Arbeiten von 
F. D. Murnaghan‘) und von B. R. Seth°), die jedoch von ganz anderen Gesichtspunkten 
ausgehen °). 


1. Die Dehnung. Nach der Formänderung werden die Koordinaten eines materiellen 
Punktes x’ zu &i. Die drei Funktionen £i stellen ein stetiges Verzerrungsfeld dar. In der 
unendlich nahen Umgebung eines Punktes x’ wird die Verschiebung d«i zu d&i: 


1) ©. B. Biezeno and H. Hencky: On the general theory of elastie stability. Proceedings of the Royal 
Academy Amsterdam. 31 (1928), p. 569, u. 32 (1929), p. 444. 

2) R.V. Southwell: On the theory of elastie stability. Phil. Trans. Roy. Soe. 1913, ser. A vol. 213, pp. 187—244. 
AR r nr refftz: Zur Theorie der Stabilität des elastischen Gleichgewichts. Z. angew. Math. Mech. Bd. 13 (1933), 

165. 

4) F. D. Murnaghan: Finite Deformations of an elastic solid. American Journal of Mathematics, 1937. 

5) B. R. Seth: Finite Strain in elastie problems. Phil. Trans. Roy. Soc. A 234 (1935), pp. 2331. 

#6) Siehe auch M. A. Biot: Theory of Elastieity with large displacements and rotations. Proceedings of the 
fifth International Congress of Applied Mechanics 1938, p. 117—122. — Thöorie de l’Elastieit& du second ordre. Ann. Soc. 
Seient. de Bruxelles LIX, ser. 1, p. 104, 1939. — Non Linear theory of Elastieity and the linearized case for a body 
under initial stress. Phil. Mag. Sec. 7, Vol. XXVII, p. 468, April 1939. 
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Zugleich besteht die Beziehung 


mit 
1 j=i 
Setzen wir 
wird 
J 
mit den Beziehungen 
i=d, wi=—al, 0. 


Das Quadrat eines Längenelementes nach der Formänderung ist 


i uv 


i 
> 


i i 


Anstatt e, haben wir e., gesetzt. 


wobei 


Der Tensor g., ist allgemein als Verschiebungstensor bekannt. Es ist jedoch besser, 
einen anderen endlichen Deformationstensor so zu bestimmen, daß er in linearer Beziehung 
zur Längenänderung ds, anstatt zum Quadrate ds? der Längenänderung steht. 

Zu diesem Zwecke wird die folgende lineare Transformation mit symmetrischen Koeffi- 
zienten eingeführt: 


wobei 
=. 


Nach Einführung dieser Transformation wird das Quadrat des Längenelementes die Form 
annehmen 


uv 

1 i 


Die Transformation (1.7) läßt bekanntlich keine Winkeländerungen zu’) und kann daher 
als Definition einer reinen Dehnung verwandt werden, wobei die endlichen Dehnungen 


durch die Komponenten &} dargestellt werden. Mit der Identität 


wobei 


wird. Wir setzen &u,= &,. 


wird die Transformation (1.4) und (1.7) den gleichen Zustand einer homogenen Dehnung dar- 
stellen. Das bedeutet zugleich, daß 

So erhält man sechs Gleichungen, in denen die sechs Dehnungskomponenten e; durch 

die Koeffizienten ei und @; der allgemeineren Transformation (1.4) dargestellt werden. Die 


?) A. E.H. Love: Mathematical theory of Elastieity, pp. 65, Cambridge 1906. 
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sechs Größen ei, die durch diese Gleichungen dargestellt sind, können als Komponenten end- 
licher Dehnung in der allgemeinen Transformation aufgefaßt werden. 

Sie sind irrationale Funktionen der Größen e! und wi: Jedoch können sie leicht be- 
stimmt werden, wenn wir uns mit einer Approximation zweiten Grades begnügen und somit 
e! und w’ als Größen erster Ordnung ansehen. 

Wir sehen aus Gl. (1.11), daß e und & sich nur um Größen zweiter Ordnung unter- 
scheiden. Wir können daher mit einem Fehler dritter Ordnung schreiben 


Dieses in Gl. (1.11) eingesetzt, ergibt 


i i 


2. Die Dehnungsenergie. Wir berechnen jetzt die Dehnungsenergie einer reinen, 
homogenen Dehnung 


J 
wobei d. 

Mit der Annahme, daß die Beanspruchung 0; und Dehnung ei Größen erster Ordnung 
sind, wird ein Ausdruck für die Dehnungsenergie hergeleitet, der bis zur dritten Ordnung 
als richtig angenommen ist. Dieses entspricht einer Näherung zweiter Ordnung der Gleich- 
gewichtsbedingungen. 

Das entsprechende homogene Spannungsfeld 0; hat die gleichen Koordinatenachsen wie 
die Dehnung. Es ist auch 


Die auf die Flächeneinheit bezogene Öberflächenkraft, die an der Oberfläche S eines 
Volumenelementes V des Körpers angreift, ist 


wobei die Größen a; die Richtungskosinusse der äußeren Flächennormalen der Oberfläche S sind. 


Ändern wir die Dehnungskomponenten um die kleinen Größen öei, so ändern sich die 
Koordinaten des Volumens um 


Die von der Oberflächenkraft geleistete Arbeit, die dieser Dehnungsänderung entspricht, 
ist gleich der Dehnungsenergie ö Wy des Volumens 


Nach (2.2) können wir schreiben 5 
Nojdfia;dS. 


Nach der Greenschen Formel kann dieser Ausdruck auch als ein Volumenintegral 


geschrieben werden: 
j i 


Da wir es mit einem homogenen Feld zu tun haben, sind die o/ alle konstant, deshalb 


ij 


Wenn man das Gleichungssystem (2.1) nach den drei Größen xi als Funktionen von &i 


auflöst, kann man die partiellen Ableitungen ade leicht erhalten. Wir lassen die Größen 


zweiter Ordnung fort und erhalten 
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i 
wo e= und D=|ö+e 


die Determinante der Transformation (2.1) ist. Diese Näherung ist gerechtfertigt, denn sie 
hat keinen Einfluß auf die Glieder der ersten drei Ordnungen im Ausdruck für öW. 


Nach Einsetzen von 
“J 


oder 
J 
in den Ausdruck (2.3), wird 


und = 
Nennen wir V’ den Volumeninhalt vor der Dehnung, so können wir schreiben 


dv 
Daher wird 


So wird denn die Änderung der Dehnungsenergie bezogen auf das ursprüngliche 
Volumen 


ij 
N N dor dei 
Der unsymmetrische Tensor 


kann daher als die Oberflächenkraft nach der Dehnung aufgefaßt werden, die an den Seiten- 
flächen eines ursprünglich kubischen Einheitselementes angreift. 

Jedoch tritt nur der symmetrische Teil dieses Tensors in dem Ausdruck von öW auf, 
weil 

Deshalb erhalten wir, wenn wir für den symmetrischen Teil den folgenden Ausdruck 
benutzen 


k 
> 
den Ausdruck: 


Der Tensor d kann als die Beanspruchung, bezogen auf die Fläche vor der Dehnung, 
und der Tensor o? als die Beanspruchung, bezogen auf die Fläche nach der Dehnung auf- 
gefaßt werden. Gl. (2.9) gibt die Beziehungen an, die zwischen den beiden Systemen von 
Spannungskomponenten bestehen. 

Das Bestehen einer Funktion der Dehnungsenergie setzt voraus, daß öW ein exaktes 
Differential ist. 

Deshalb müssen die Spannungs-Dehnungs-Beziehungen 


.&..) oder 
den folgenden achtzehn Bedingungen genügen: 


dr 
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oder 

1 ‚da ; _ (of + 


Hieraus ist klar ersichtlich, daß die Beziehungen für das Spannungssystem r nicht die- 
selben sind wie für das System o. 


3. Die Gleichgewichtsbedingungen. Die Änderung der Dehnungsenergie, bezogen auf den . 
Volumeninhalt V für den Fall der niehthomogenen Dehnung ist 


uw 


Die Größen €; müssen alle Glieder erster und zweiter Ordnung enthalten, so, wie sie 
vorher berechnet worden sind: 


i i 
1 1 
= Eur Eur +z (wi, etw ei) + 3 > 


Da in unserem Falle nur cartesische Tensoren benutzt sind, können für jeden Tensor A, 
auch die Bezeichnungen A”“ oder A». benutzt werden. 
Durch Variation ergeben sich die folgenden Ausdrücke: 


2 2 u 


Zudem kann unter Benutzung der Beziehungen (1.3) und durch richtigen Gebrauch der 
Summationszeichen ag werden, daß 


Wir benutzen die Eigenschaft, daß Ps 2 öui und integrieren partiell. Dann er- 
halten wir aus dem Integral (3.1) die Beziehung 


||| 


v uv 


und worin a, die ee der äußeren Normalen der Oberfläche von V sind. 


Die virtuelle Arbeit der inneren Kräfte ist gleich — ö5W und die virtuelle Arbeit der 
äußeren Kräfte ist gleich: 


° 


ist. 


wo 


2.2... 
wobei Fi die Oberflächenkraft ist, bezogen auf die ursprüngliche Fläche, und Xi die Massen- 


kraft darstellt. Wenn wir bedenken, daß die gesamte virtuelle Arbeit für alle Werte von öwi 
verschwindet, so erhalten wir 


und die Gleichgewichtsbedingungen nehmen die Form 


v uv 
ri 


94 Biot, Elastizitätstheorie zweiter Ordnung mit Anwendungen AA 


an, und die Randbedingungen werden durch 


dargestellt. 
Diese Gleichungen können auch durch die Spannungen o“” dargestellt werden. Da wir 
uns auf Glieder erster und zweiter Ordnung beschränkt haben, können wir für &, die Größe 


e, einsetzen. 


u 


Diese Ausdrücke werden in die Gl. (3.7) und (3.8) eingesetzt. Es genügt, wenn wir für 
unsere Näherung r’i nur in die Glieder erster Ordnung einsetzen und in den übrigen Gliedern 
an Stelle von ri die Größen 0”! nehmen. 

Auf diese Weise erhalten wir 


und die Randbedingungen 


v uv 
Fi= N (l+e)oria, + N — 
4. Die lineare Elastizitätstheorie eines vorbeanspruchten Körpers. Als Anfangszustand 
nehmen wir an, daß eine Vorbeanspruchung S“” im Körper existiert, der die Koordinaten «'i 
entsprechen. Die Gleichgewichtsbedingung dieses Anfangsspannungsfeldes lautet: 


worin o Xi die anfängliche Einheitsmassenkraft bedeutet. 

Nehmen wir jetzt an, daß der Körper einer kleinen Dehnung unterzogen wird. Für 
diese Dehnungen bauen wir eine lineare Elastizitätstlieorie auf. 

Die Koordinaten xi eines materiellen Punktes im Körper werden durch die Dehnung zu 
xzi+ ui. Die Koordinatenänderungen wi und der dazugehörige Tensor e;, »; werden als 
Näherungen erster Ordnung angesehen. Die Dehnungskomponenten e; sind wie in Gl. (1.13) 
definiert. Für die Komponenten der Spannung nach der Dehnung können nun entweder die 
tatsächliche Spannung S""+ oder die Spannung bezogen auf die 
Fläche vor der Dehnung gewählt werden. Die Spannungsänderungen #“” und s“’ sind auch 
als klein und als Näherungen erster Ordnung anzusehen. 

Die einer linearen Theorie entsprechende Dehnungsenergie muß nun lineare und qua- 
dratische Glieder enthalten. Wir behandeln somit die Änderung der Dehnungsenergie, öW, 
die einer Änderung der Dehnung von der Größenordnung ö &., entspricht. Man erkennt, daß 
es möglich ist, einen Ausdruck für öW zu erhalten, der alle Glieder erster und zweiter 
Ordnung enthält, wenn man 0“” = S“" + s“” in die Gl. (2.10) einsetzt. Wenn man alle Glieder 
dritter Ordnung vernachlässigt, erhält man 


uw uv 
worin 
ist. 


Die obigen Gleichungen drücken die Beziehungen zwischen den beiden Systemen der 
Spannungskomponenten s“” und #“” aus: Die Spannungs-Dehnungsbeziehungen sind als linear 
angenommen. Die Form dieser Beziehungen hängt natürlich von der Wahl des Spannungs- 
systemes ab, d.h. ob wir das System s“” oder {“” benutzen. 


“r 
oder 


“r 


AUM 
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“r 
Die Summierung I’ wird über alle sechs Kombinationen von x und r ausgedehnt. Da- 
mit öW ein totales Differential sein kann, müssen gewisse Beziehungen zwischen den 
Koeffizienten B und € bestehen. Es bestehen die Gleichungen 


(4.6), 
somit 

ET 

Bir + +S ug, + 59’; 
(4.8). 

dE 1 


Die C-Koeffizienten stellen eine symmetrische Matrix dar, während dasselbe für die B- 
Koeffizienten im allgemeinen nicht zutrifft. 


uv 
Der Ausdruck öW’ = N’ t“"ö eu, stellt das Differential einer quadratischen Form W’ in 


der Veränderlichen e., dar. i : Hr Se und von den Eigenschaften quadratischer 
Formen 
w 


So erhalten wir einen Ausdruck für die potentielle Energie 


uv uv 
1 


der alle Glieder erster und zweiter Ordnung enthält. 
Weiterhin kann an Stelle von eu, auch e., eingesetzt werden, zum mindesten in der 


uv 
ersten Gruppe der Glieder I’#“" &u,, denn W wird dadurch nur in der dritten Ordnung beein- 


flußt. Somit 
1 uv 
W=z ) eur + ) SUr gu» (4.11), 


wobei die Beziehung zwischen Spannung und Dehnung durch 


kr 
gegeben ist. 


Zugleich können auch die Spannungen s“” aus (4.3) eingeführt werden und man erhält 


k 


Nach Einsetzen dieser Ausdrücke inW erhalten wir 
S’k + ei. S"k) . (4.14) 
mit den Spannungs-Dehnungsbeziehungen 
kr 


Es muß nun klar im Auge behalten werden, daß der Wert &.,, der hier benutzt wird, 
durch die Gl. (1.13) dargestellt ist: 


i 
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Wenn wir nun diese Ausdrücke der Dehnungsenergie W benutzen und die Variations- 
methode der virtuellen Arbeit anwenden, wie es im nichtlinearen Falle geschehen ist, sowie 
die Gleichgewichtsbedingung (4.1) des Vorspannungsfeldes S“” berücksichtigen, so erhalten 
wir die folgenden Gleichungen für das Gleichgewicht 


v uv 
pi 1 


worin 4 Xi die Änderung erster Ordnung der Massenkraft ist. 
Wir erhalten für die N 


worin 4 Fi die Änderung der Oberflächenkraft, bezogen auf die Einheit der ursprünglichen 
Oberfläche, und a, der Kosinus der Richtungswinkel der nach außen gerichteten Normalen 
bedeuten. 

Wir haben somit die Gleichungen für die Spannungsänderung #“”, bezogen auf die 
ursprünglichen Flächen. Wir können auch die tatsächliche Spannungsänderung s“” einführen, 
entweder indem wir mit dem Ausdruck für W (4.14) beginnen, oder indem wir die Beziehungen 
(4.13) in die Gl. (4.16) und (4.17) einsetzen. 


Wir erhalten so die Gleichgewichtsbedingungen 


und die Randbedingungen: 


v v uv 


Die Gl. (4.16), (4.17), (4.18), (4,19) enthalten alle Glieder erster Ordnung, die in einer Theorie 
der Elastizität eines Körpers unter Vorspannung enthalten sind. 

Die Gl. (4.18) können auch in einer anderen Form geschrieben werden, wenn wir die 
Gleichgewichtsbedingungen (4.1) berücksichtigen, die sich auf die Vorspannungen beziehen, 
sowie die Identitäten: 


u 
de, de 


u 
worin e= Ne, 


Die Gl. u. haben die Form 


i 


Die Glieder in Gl. (4.20) lassen eine interessante Deutung zu. Nehmen wir zum Bei- 
spiel den Fall einer zwei-dimensionalen Dehnung und setzen 


worin 
_du 


Mit diesen Bezeichnungen erhalten die Gleichgewichtsbedingungen (4.18) die Form 
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1 
+ FedX+ooY 

REM) dm 

+(S Ss —28 


und eine ähnliche Gleichung für die Y-Richtung. 

Die zusätzlichen Glieder in der zweiten und dritten Reihe können als Folge der Form- 
änderung eines infinitesimalen Elementes aufgefaßt werden. Die Glieder in der zweiten Reihe 
stellen den Einfluß der Krümmung dar. Die Glieder in der dritten Reihe haben ihren Ursprung 
in einer Art Auftrieb eines Elementes durch seine Formänderung im ursprünglichen Dehnungs- 
feld. Man erkennt somit, daß die Glieder, die die Krümmung darstellen, nur von ursprüng- 
lichen größten Querspannungen abhängen und für den Fall eines anfänglichen gleichmäßigen 
Druckes verschwinden. Die Glieder, die mit dem Auftrieb zu tun haben, verschwinden, wenn 
das ursprüngliche Spannungsfeld homogen ist. In diesem Falle, oder wenn wir den anfäng- 
lichen Spannungsgradienten und die Massenkraft vernachlässigen, erhalten wir die Gleichungen 
—S dy 28 = =0, de +(8 Ss dy —=0 (4.22). 

Diese letzte Art von Gleichungen muß eine Ben und fundamentale Rolle in Knickungs- 
aufgaben spielen. Man erkennt, daß überhaupt die Möglichkeit einer Knickung im wesent- 
lichen von der Existenz einer Querspannung im ursprünglichen Spannungsfeld abhängt. Diese 
Folgerungen sind auch im drei-dimensionalen Falle richtig. 


5. Anwendungen. Die Dehnungskomponenten zweiter Ordnung in der 


Plattentheorie. Die Verschiebungen in Biegung und Dehnung werden angenommen 
wie folgt: 


wenn). 
Die Ebene z2=0 wird als Mittelebene vor der Biegung angenommen. Durch Anwendung der 
Formel (1.13) erhält man die drei wichtigsten Komponenten des Dehnungstensors: 
(3= 
Isotrope Beziehungen zwischen Spannung und Dehnung mit ENT 
tischen Gliedern. 
Wegen der Isotropie fallen die Hauptachsen für Spannung und Dehnung zusammen. 
Wir wählen die rechtwinkligen Achsen 1, 2, 3, parallel mit den Hauptrichtungen und be- 
trachten jetzt die Spannungskomponente r,, als Funktion der Dehnungskomponenten £,,, #ge, Ess; 
die Glieder erster und zweiter Ordnung enthalten. Ein Vertauschen der Achsen 2 und 3 hat 
wegen der Isotropie keinen Einfluß auf r,,, deshalb 


In ähnlicher Weise erhalten wir: 


worin €&=&,,+#3-+ £;, die kubische Dilatation bedeutet. Unter der Annahme eines elastischen 
Potentiales müssen die Beziehungen (2.11) befriedigt sein. 


Leitet man noch die Beziehung 2B=( ab, werden die Beziehungen zwischen Spannung 


und Dehnung durch fünf Elastizitätskonstante dargestellt. Das Bestehen von fünf physi- 
kalischen Konstanten war auch schon von F. D. Murnaghan festgestellt worden, der aber 


(5.1). 


7 
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auf anderem Wege zu diesem Schluß gelangt ist (siehe oben). Die physikalische Bedeutung 
seiner Konstanten ist jedoch eine andere, als in unserem Falle. 
Die Spannungs-Dehnungsbeziehungen können in der äquivalenten Form dargestellt 
werden: 


Im Falle, daß wir es nicht mit den Hauptrichtungen zu tun haben, erhalten wir in Carte- 
sischen Koordinaten 


aß @ 
Le? + N N 803809) gur+ (u+2Ne)eu+ M (53), 


worin Qu»=1, wenn u=r und Qu,=0, wenn ur. Diese fünf Elastizitätskonstanten A, G, 
L, M, N könnten bequem zu einer genaueren Beschreibung sehr elastischer Stoffe benutzt 
werden. 

Erhöhung der Drehungssteifigkeit eines prismatischen Stabes durch 
axiale Spannung. 

Der prismatische Stab wird unter gleichmäßiger axialer Spannung S,, = S angenommen; 
die z-Achse ist parallel mit der Stabachse. Als zusätzliche Spannungskomponenten der Drehung 
werden die s«, angenommen. Die Gleichgewichtsbedingungen sind (4.20) 


0815 


dx + + dz dz . . . . . . . . (5.4). 


Die Randbedingungen (4.19) werden 
AF, +s,0,+Soya,, 
AF,„= + 830+ 8, a, — So; a,, tar; 
AF,=s,0, ,+s,0,+S(e, + e,) — Se,a, —Se,a, 


Als angenähertes Spannungs-Dehnungsgesetz wird das Hookesche Gesetz für ein iso- 
tropes Material angenommen: 


5.6). 


Es kann leicht gezeigt werden, daß die klassische Lösung der Torsionsaufgabe von Saint 
Venant, die sich auf einen Stab ohne Vorspannung bezieht, auch die obigen Gleichungen 
befriedigt. Wenn © den Drehungswinkel bezeichnet, dann bedeutet s,, und s,, die Verteilung 
der Querbeanspruchung über die Schnittfläche und ist durch die folgenden Gleichungen dar- 
gestellt: 


085, , OS _ 


wobei die Randbedingungen aussagen, daß die Querspannung tangential zu der Schnittfläche 
sein muß. Wenn wir nun das an der Schnittfläche angreifende Drehmoment berechnen, er- 
halten wir nicht dasselbe Resultat, das die Saint Venantsche Theorie ergeben würde. Das 
hat seinen Grund in der Tatsache, daß die an der Schnittfläche angreifenden Kräfte, die durch 
die Randbedingungen (5.5) dargestellt sind, die Form haben: 


AF,=s8,+Sw,, =, + 
Die Werte von », und &, werden aus der Saint Venantschen Lösung abgeleitet, und wir 


erhalten: 


Das an der Schnittfläche angreifende Gesamtdrehmoment ist: 


T=(1--35)\ \ 


XUM 
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oder 
T=(1-55) Tr . 


In diesem Ausdruck bedeutet Ty;y das aus der Saint Venantschen Theorie abgeleitete 
Drehmoment und I, das polare Trägheitsmoment der Schnittfläche, bezogen auf den Schwer- 
punkt der Fläche. 

Der Nullpunkt des Koordinatensystemes x, y muß hier mit dem Schwerpunkte zusammen- 
fallen, damit das aus diesen Gleichungen erhaltene Kraftsystem ein reines Drehmoment ist. 


Das Glied . kann im Vergleich mit „Eins“ vernachlässigt werden, weil es eine Näherung 


gleicher Ordnung ist, wie die Näherung durch die Einführung des Hookeschen Gesetzes für 
ein isotropes Material (5.6). 
Deshalb ist der Wert des Drehmomentes 


Das Korrekturglied für die anfängliche axiale Vorspannung erhält nur dann Bedeutung, wenn 
die Drehsteifigkeit der Saint Venantschen Theorie klein ist im Vergleich zu S I;. Bei- 


Bild 1. 


spielsweise erhalten wir im Falle eines dünnwandigen Zylinders, der parallel seiner Achse 
aufgeschnitten ist, 


2 


wobei „r“ den Radius und „c“ die Wandstärke bedeuten. Wählen wir zum Beispiel r/c = 10 
und S/@ =1/1000, so erhalten wir durch den Einfluß der axialen Spannung einen Fehler des 
Saint Venantschen Drehmomentes von 30°s. 118 


Note on the bending of thin uniformly loaded Plates bounded 
by Cardioids, Lemniscates and certain other Quartic Curves. 
By Bibhutibhusan Sen, Krishnagar College, Bengal. 


Zusammenfassung. Die Formeln für die Durchbiegung dünner, gleichmäßig belasteter 
Platten, die von Kardioiden, von Lemniskaten oder von Kurven begrenzt sind, die zur Ellipse 
bzw. Hyperbel in bezug auf den Mittelpunkt oder zur Ellipse in bezug auf den Brennpunkt 
als Pol (elliptie Limacon) invers sind, werden abgeleitet. Für den Fall, daß die Ränder ein- 
geklemmt sind, kann die Lösung in geschlossener Form gegeben werden. 


Introduction. The normal displacement w of a thin plate bent by transverse forces, 
satisfies the differential equation ') 
%w 
where Z=the constant load per unit of area, and 
D=the flexural rigidity. 


For a clamped edge, we have on the boundary 


dv being an element of outward drawn normal. 


1) A.E.H. Love: The Mathematical Theory of Elastieity. 3rd. edition (19%), p. 492. 
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The exact solution of the equation (A) satisfying the boundary conditions (B) is known 
only for some thin plates of simple forms. The object of this paper is to consider the same 
problem when the boundaries are Cardioids, Lemniscates, Inverses of ellipses and elliptie 
Limacons. The results are interesting in so much as they are expressed in finite number 
of terms. 


1. Solution for a Plate bounded by a Cardioid. By using the transformation 


1 
we obtain 
=r7?70087, sinz, 
and 


It is evident that constant and 7 constant give two families of orthogonal Cardioids. 


We shall take as the boundary of the plate 
f=e= - (2a): r=all+c9). ........ (13. 


The above transformation enables us to write the equation (A) in the form 
E-a 
where 
15 
8 (1.5). 
ä The particular a of this equation is given by 
Combining this expression with some suitable terms which make 
0° 2 
we obtain 
(1.8), 
where P,, P,, P,, P, are constants. 
From the boundary conditions (B) we get 
when !=a 
These relations give 
When these values of the constants are substituted in (1.8), we have 
Expressed in polar coordinates, this result stands as 2 
5 
(1.11). 
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It is apparent that the above expression for the normal displacement is symmetrical about 
the line 9=0. The value of the defleetion along this line is given by 


= ae 1 -— 2... 
where o is the distance measured along the line of symmetry from the origin. 
On putting 
e 


the above expression becomes 


By simple caleulation we find that the maximum displacement on the line of symmetry oceurs 
at a distance 


where its approximate value is 
8... (1.16). 


2. Solution for a Plate bounded by a Lemniscate of Bernoulli. Let us consider the system 
of orthogonal coordinates (£,n) defined by 


where (2.1). 
and z=r-+iy 
In cartesian eoordinates we obtain the relations 


so that £= constant gives a family of quartic curves?) with real foei at z=+c and „= con- 
stant gives the orthogonal family of curves which are rectangular hyperbolas, each of which 
passes through the foci of the quartic curves. From (2.2) it is evident that the curve E=0 
gives a Lemniscate of Bernoulli of which the polar equation is 


At infinity while at either focus x. We shall consider 
the loop of the Lemniscate of Bernoulli given by the curve Z=0 £-0 
which lies in the region 
s \m 
2 2 6 
Since in the present case 


the equation (A) can be expressed in terms of £ and n in the form 


(26) 

The particular integral of the above equation will be 


- Finding out a few terms which make the left hand side of (2.5) equal to zero, we can write 
for the purpose of our problem 


P,e°cosn (cosh2& + 
where P,,P, and P, are constants. 


2) Basset: Treatise on Hydrodynamices, Vol. 1 (1888), p. 106. 
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For a clamped edge we must have on the boundary 


These conditions are satisfied if 


Hence the normal displacement in this case is obtained as 


. . . .@.M). 
At the focus where > — x 
3. Solution for a Plate bounded by an Inverse of an Ellipse. If we use the transformation°) 
we get 
_.2ecoshncos$ .2esinhnsin& 
 cosh2n +cos2£’ 9 + cos2E > 
(3.2) 
and 
2.c? (cosh 27 — cos 2£) 38 
(cosh2n + cos2 em 


It is evident that is the inverse of an ellipse and £=a is an inverse of a confocal 
hyperbola which becomes an arc of a Lemniscate of Bernoulli when e=T. Here £ may 


have any real value positive or negative but as the values of x and y are periodie in £&, it 
is only necessary to consider values of £1ying between 0 and 27. Values of n are restrieted 
to all positive real values between zero and infinity. At the origin „>x and when n=0 
we have two double lines extending from the points <=+c to infinity in the positive and 
negative directions respectively. 


In the present case the equation (A) anllanee to 
else _4eZ 


—=64p (say) (3.4). 


cosh 27 —cos2£ D 

For clamped edge the boundary conditions are 

N 
7 for an Inverse of an ellipse. For an Inverse of a hyperbola 
we shall have 


As the solution of the equation (3.4) we can take‘) 


_P+P, (cosh2n cos22+1)+P,sinh2n 
(cosh 27 + cos 2 £)? 


where P,, P,, P,, P, are constants. 


P,sinh2n+P, 
cosh2n+cos2E 


+ (8.7), 


3) Basset, p. 102. 
4) For the Inverse of a hyperbola we can take the solution of (3.4) as 


p + 0Qı (eosh 2 ncos2& +1) + Q,sin25+Q 


(eosh 27 + c08 25)? 


+ 


8.7), 


where Q,, Q2, Qı are constants. 
Provided £ be not equal to z/4 or a multiple of 7/4, the boundary conditions (3.6) will give 


Qı=peosec?2«, Qz=—2peosee 2a, 
Qs=2peosec4a, Qu =—2p eosee 2a cosee 4a 
Special solutions can also be obtained when #=7/4 or a multiple of n/4. 
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The boundary conditions (3.5) will be satisfied if 


P,=2p cosech 4 ß, P,=2p cosech 2 ß eosech 4 ß si 
At the centre where „>, we have 


4. Solution for a Plate bounded by an elliptice Limacon. By using the transformation °) 

e 
we find that 7 = ß (constant) represents an inverse of an ellipse with respect to the focus as 
the pole, that is, an elliptic Limacon. The values of 7 are restrieted to all real positive 


values between zero and infinity while £ may have any real value positive or negative. 


From the above transformation we get 


These relations enable us to write the equation (A) in the form 
(cosh 7 + cos [ Z_ 
The solution of this equation suitable for the present problem can J 
be obtained as 
p+P, (cosh? n cos? &—1)+ P, (cosh cos&+1) 7B, 
w= 
(cosh n + cos £&)* re | 
(cosh n + cos £)? (cosh n + cos &)? 
P,,P,,P,, P,,P,, P, being constants. 
On the boundary n=f, we have Bild 4. 
These conditions are satisfied if 
P,=p cosech* ß, P,=2P, cosh? ß, P,=— coth B(4sinh’ ?+3), 
> (4.6). 
P,=5'coth ß, P,=-5tanh2ß, P,=— P, cosh? sech 2 8 
At the origin where >», we have 
cosech* . . . 2.2.2.7. 


Summary: In this note the bending of thin uniformly loaded plates bounded by Cardioids, 
Lemniscates and certain other quartice curves is discussed. Assuming that the boundary is 
elamped, the solution in each case is obtained in finite number of terms. ' 124 


5) Basset, p. 106. 
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Zur numerischen Integration gewöhnlicher Differential- 
gleichungen zweiter und höherer Ordnung”). 


Untersuchungen zu den Verfahren von Blaeß und Runge-Kutta-Nyström. 
Von Rudolf Zurmühl in Darmstadt. 
(Aus dem Institut für Praktische Mathematik der Technischer Hochschule Darmstadt, Prof. Dr. A. Walther.) 


1. Ziel der Arbeit. Für die numerische Integration von Differentialgleichungen ') zweiter 
und höherer Ordnung gibt es, trotz einer großen Anzahl allgemeiner Lösungsmethoden?), nur 
wenige Verfahren, welche den Näherungsaufbau in unmittelbarer Form, ohne ein Auf- 
spalten in ein System von Gleichungen 1. Ordnung, vornehmen. Wenn sich trotz der 
Wichtigkeit namentlich der Differentialgleichungen 2. Ordnung für die technischen Wissen- 
schaften die Behandlung durch numerische Integration in der Ingenieur-Mathematik noch 
nicht so wie erwünscht und möglich vorfindet, so liegt das zu einem Teil in dem Mangel an 
unmittelbaren und damit auch handlichen Lösungswegen begründet. 


Die Verwendungsmöglichkeit numerischer Integrationsverfahren für Ingenieuraufgaben 
ist noch bei weitem nicht ausgeschöpft. Das technische Schrifttum hat eine Fülle von teils 
recht verwickelten Sondermethoden in Fällen hervorgebracht, in denen mit einem einzigen 
brauchbaren Verfahren unmittelbarer Integration der zugrunde liegenden Differentialgleichung 
eine bessere und meist auch raschere Lösung erzielt wird?). 


Zu den direkten Näherungsmethoden für Differentialgleichungen höherer Ordnung gehört 
das Verfahren von Blaeß*), welches insbesondere für Differentialgleichungen 2. Ordnung 
einen sowohl numerisch wie graphisch gangbaren Lösungsweg angibt. Die in ihm entwickelte 
Rechnungsform ist der Eigenart der Differentialgleichung 2. Ordnung formal wie anschaulich 
besonders angepaßt. Sie kann daher über die Grenzen des Verfahrens selbst hinaus von 
Nutzen sein. 


Die vorliegende Arbeit setzt sich, nach einer kurzen Einführung in das genannte Ver- 
fahren, zwei Aufgaben: 
1. Untersuchung über gewisse Fragen der Genauigkeit der Blaeßschen Näherung. 


2. Unmittelbare Behandlung von Differentialgleichungen höherer Ordnung nach dem 
Prinzip des Runge-Kutta-Verfahrens. Für Differentialgleichungen 2. Ordnung 
ergeben sich dabei bekannte Formeln von Nyström°) in neuer Fassung. Für 
Gleichungen 3. Ordnung lassen sich neue Formeln herleiten. 


I. Das Verfahren von Blaeß. 


2. Prinzip der Näherung. Das Verfahren fußt auf dem Taylorschen Satz und gehört zu 
den Lösungsmethoden, welche eine Näherung durch möglichst weitgehendes Angleichen an 
eine Taylorentwicklung aufzubauen suchen °). 

Es ist bekannt, daß zu einer gewöhnlichen Differentialgleichung in ft und y aus ge- 
gebenen Anfangswerten #,, %,, %,, ... ein Näherungswert y, an einer um einen Schritt dt=h 
entfernten Nachbarstelle , =t,+h nach der Taylorentwicklung berechnet werden kann zu 


h* 


*), D 87 (gekürzt). 

»)H. von Sanden: Praktische Analysis, Berlin 1923, S. 162 bis 190. — C. Runge, HA. König: Vorlesungen 
über numerisches Rechnen, Berlin 1924, S. 286 bis 323. — Fr. A. Willers: Methoden der praktischen Analysis, Berlin 
1928, S. 295 bis 339. 

2)C. Runge, Fr. A. Willers: Numerische und graphische Integration gewöhnlicher und partieller Differen- 
tialgleichungen. Enz. d. Math. Wiss. Bd. II, C 2, S. 147 bis 150. — A. A. Bennet, W.E. Milne: Numerical Inte- 
gration of Differential Equations. Bull. of the National Research Couneil, Nr. 92. Washington 1933, S. 71 bis 78. 

3) Wie weit diese Behauptung zutrifft, zeigt ein von V. Blaeß angegebenes Verfahren zur Lösung von Rand- 
und Eigenwertaufgaben durch numerische Integration, wobei der Eigenwert als unbestimmter Parameter durch die 
Rechnung mitgenommen wird: V. Blaeß: Zur Berechnung der Schwingungen stabförmiger Körper mit veränder- 
lichem Querschnitt. Ing.-Arch. Bd. 9 (1938), S. 428 bis 435. 

4) V.BlaeBß: Zur angenäherten Lösung gewöhnlicher Differentialgleichungen. Z. VDI Bd. 81 (1937), S. 587 bis 596. 

5) E.J. Nyström: Über die numerische Integration von Differentialgleichungen. Acta Soc. Sei. fennieae Bd. 50 
(1925), Nr. 13. 

6) Vgl. Runge-Willers:a. a. O.?), S. 147 bis 150, und Willers: a. a. O.!), S. 305 bis 314. 
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Für eine Differentialgleichung 2. Ordnung 


sind die beiden ersten Glieder durch die Anfangswerte y, und ,=v, an der Anfangsstelle 
t=t, gegeben. Das nächste Glied ij,h’/2 läßt sich durch Einsetzen der gegebenen Anfangs- 
werte in 5, =f(t,9.,%) unmittelbar aus der Differentialgleichung berechnen. 
Alle folgenden Glieder mit ihren höheren Ableitungen müssen durch Weiterdifferen- 
zieren der Differentialgleichung (1) gewonnen werden. Da dieses Weiterdifferen- 
zieren, weil es unter Berücksichtigung aller in der Differentialgleichung auftretenden Ver- 
änderlichen erfolgen muß, in Gl. (1) also hinsichtlich £, y und y=», in den meisten Fällen 
sehr mühsam, praktisch sogar bald unmöglich wird, so ist dieses Verfahren der Taylor- 
entwicklung als selbständiges Näherungsverfahren kaum in Gebrauch ’). 

Das Verfahren von Blaeß bedient sich des Taylorsatzes in unmittelbarer Form 
gerade bis zu dem Gliede, welches aus der Differentialgleichung durch einfaches Ein- 
setzen, ohne Weiterdifferenzieren zu berechnen ist. Für eine Differentialgleichung 
2. Ordnung also wird die Taylorreihe abgebrochen nach dem Glied mit 7. Die so erzielten 
Näherungen, die wir mit 7 und 7=®# bezeichnen wollen, gelten als Rohwerte. Sie werden 
verbessert durch einen gesonderten Fehlerabgleich in Form von Korrekturen; er führt 
die Einbeziehung weiterer Taylorglieder herbei, ohne i. a. ein Weiterdifferenzieren der 
Differentialgleichung zu benötigen. 

Für eine Differentialgleichung 2. Ordnung lautet die Berechnung der Rohwerte %, und 
®, nach der verkürzten Taylorformel, geschrieben in der mit 4 dimensionsgleichen Tan- 
gentenerhebung »h und mit der gleichfalls auf y-Dimension gebrachten „Abwölbung“ 


folgendermaßen: 


Die numerische Rechnung folgt dem Schema Tafel 1°), welches, außer der Berechnung der 
aj;-Werte durch Einsetzen der Werte t;, %, ®; in die Differentialgleichung, nur einfache Ad- 
ditionen erfordert. Die Schrittrechnung des Schemas wird nach einer Anzahl von Roh- 
schritten, zweckmäßigerweise fünf, durch die (Grob-)Korrekturrechnung (siehe 3.) unter- 
brochen, welche durch Hinzufügen gewisser Verbesserungen V und V’ zu den Rohwerten 
%, und ®,h bessere Näherungswerte 4, und v,h liefert. Die Rechnung erfährt so eine Unter- 
teilung in Gruppen zu je fünf Rohschritten. Anfangswerte jeder neuen Gruppe sind die 
korrigierten Endwerte der alten. 


Tafel 1. Rechnungsgang nach Blaeß. 


Nr. t Yi vih a; 

0 0 Yo %)-h2/2 

1 h | a, =f(h, 

2 2h | ,—/(2h, 

3 3h =... 

4 4h =... | 

Korr. V V' 

5 5h 

6 6h 


Anschaulich bedeutet die Näherung nach den Schrittgleichungen (3) den Ersatz der 
strengen Lösungskurve y(f) mit ihrer gemäß der Differentialgleichung stetig veränder- 
lichen Wölbung %=f(t,y,v) durch einen „Parabelzug“ %(f) von stückweise fester, 
jeweils zum Schrittanfang nach ij; =f(t;, errechneter Wölbung. Der auf 4-Dimension 
gebrachte Ausdruck ; =; h?j2 stellt dabei die Abwölbung des Parabelbogens %(f) am 
Schrittende von der Anfangstangente dar. Die im Zahlenschema verfolgbare Einwirkung der 
mit #,%,% laufend sich ändernden a-Werte auf die Steigungen v gewährt einen unmittelbaren 
Einblick in den Mechanismus der Differentialgleichung 2. Ordnung, weiche den Verlauf 
der Lösungskurve durch die mit £,9,v veränderliche Wölbung steuert. 


Siehe Bennet-Milne: a. a. O.?2), S. 71 bis 72 oder Willers: a. a. O.!), S. 305 bis 306. 
8) Es beginnt mit to =0. 
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Den geometrischen Aufbau 
der Näherungswerte %, und #,h 
aus den Anfangswerten 4, und v, h 
und aus der zum Schrittanfang 
nach Gl]. (1) errechneten Abwöl- 
bung a, veranschaulichen die Bilder 


1 und 2. 


Abgesehen von den Besonder- . 

| heiten des Fehlerabgleichs, auf den 
ji % % noch einzugehen ist, sind für die 

x praktische Handhabung des Ver- 

fahrens wesentlich: 
Bild 1 und 2. Aufbau der Näherung nach Blaeß. zum Näherungsaufbau; 


2. die Durchführung der Zahlenrechnung in einheitlicher y-Dimension. 


Beide Merkmale führen Anschaulichkeit und Einfachheit der Zahlenrechnung herbei. 
Von beiden läßt sich auch bei anderen Lösungswegen vorteilhaft Gebrauch machen. 


3. Der Fehlerabgleich. Die Verbesserungen V und V’, die den rohen Gruppen- 
endwerten %, und ®,h zuzuschlagen sind, um eine nachträgliche Einbeziehung 
höherer Taylorglieder in die Näherung ohne ein Weiterdifferenzieren der 
Differentialgleichung zu ermöglichen, werden nach Blaeß als Jineare Verbin- 
dungen aus den im Verlaufe der Schrittgruppe auftretenden veränder- 
lichen a;-Werten berechnet. Sie lauten: 


45 100 145 5 16 


Die damit erhaltenen Näherungen am Gruppenende, zugleich die Anfangswerte der neuen 
Fünfergruppe, sind: 


Den Werten V und V’ in Gl. (4) liegt, wie erforderlich, die Veränderlichkeit der Werte %; 
(anschaulich der Abwölbungen «;) innerhalb des Gruppenschrittes 5h zugrunde. Beide Aus- 
drücke V und V’ werden Null für y=konst, wenn also höhere Taylorglieder fehlen. Auf die 
Herleitung der Formeln und die Auswahl der in ihnen herangezogenen a-Werte wird im 
II. Abschnitt eingegangen. 

Außer der soeben besprochenen (Grob-) Korrektur, welche für. mittlere Genauigkeiten, 
also namentlich bei technischen Anwendungen, fast immer ausreicht, ist von Blaeß noch 
eine weitergehende „Feinkorrektur“ entwickelt worden®), mit der sich im Bedarfsfalle 
die Genauigkeit wesentlich steigern läßt. Die Feinkorrektur verwendet die ersten partiellen 
Ableitungen f, und der Differentialgleichung =f(t, y,v). Sie wird jedoch in der vor- 
liegenden Arbeit nicht berücksichtigt. 


II. Berechnungen zum Fehlerabgleich nach Blaeß. 


4. Ziel der Rechnung. Zum Zwecke der Rechnungsvereinfachung und kurzen Dar- 
stellung sind a. a. O.*) die zur Herleitung der Korrektur erforderlichen Differentiationen und 
Taylorentwicklungen allein hinsichtlich t angegeben worden'°).. Hier wird eine Durch- 
rechnung unter Beachtung aller drei Veränderlichen #£, y, v gegeben. 

Zu vergleichen sind die rohen Gruppenendwerte %, und ©, h mit den Werten %, und ©, h einer 
Taylorentwicklung vom Gruppenanfang her über einen Gruppenschritt Der Unter- 
schied von Taylorwert gegen Rohwert, also , —%, bzw. %,h—ü,h, den wir den „Fehler“ 
nennen wollen, soll durch Verbesserungen V und V’ (Korrekturen), welche den Rohwerten zuzu- 
schlagen sind, soweit wie möglich zum Abgleich gebracht werden. 

Die Verbesserungen V und V’ werden dabei a's lineare Verbindungen aus den a-Werten 


» der Gruppenrechnung aufgebaut, also nach einem Ansatz: 


(6). 


V=4%+4 +49, +4; +49, +40 
V’=Ba+ 


®) Bisher nicht veröffentlicht. Vgl. aber Blaeß, a.a.0.%), S. 589 links, 18. bis 15. Zeile von unten. 
10) Blaeß, a. a. O.*%),.S. 589, links: „An Stelle von nur drei Gliedern ..... 
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Darin sind die noch unbestimmten Beiwerte A und B so zu wählen, daß zwischen „Fehlern“ %, — %; 
bzw. %,h— ©, h und Verbesserungen V bzw. V’ eine möglichst weitgehende Übereinstimmung hinsichtlich 
der Taylorglieder erreicht wird. 


Eine Taylorentwicklung über den Gruppenschritt 4£=5h liefert für 4, und v,h die Taylor- 
werte: 


Die rohen Gruppenendwerte % und ?,Ah gehen aus den fünfmal hintereinander angewandten 
Schrittgleichungen (3) hervor zu: 


Ein Vergleich der Werte (7) und (8), also das Aufstellen der Fehler, sowie das Nachbilden der Fehler 
durch die Verbesserungen V und V’ nach Ansatz (6) setzt zunächst voraus: 


1. Umschreiben der höheren Ableitungen %,y!Y, yV,... auf Differentialausdrücke 
derzweiten Ableitung ä#=f(t,y,v) =} durch vollständiges Differenzieren der Differential- 
gleichung (1) nach t; 


2. Entwicklung der Werte a, bis a,, also /, bis /, in eine Taylorreihe aus dem Ausgangs- 
wert /—=/ und dessen Differentialausdrücken. 


5. Durchführung. Die höheren Ableitungen /, yIV, yY, ... entstehen durch mehrfaches 
vollständiges Differenzieren der Differentialgleichung #—=/(t,y,v) nach {. Man führt hierzu zweck- 
mäßig einen Operator D ein, welcher, angewandt auf eine beliebige Funktion 9 (£,y,v), den Ausdruck 


liefert mit und 


Außer den gebräuchlichen Kürzungen 


mit symbolisch verstandener Klammerpotenz werden gewisse Ausdrücke, welche sowohl bei den 
höheren Ableitungen als auch bei der Taylorentwicklung der a;-Werte in gleicher Form auftreten, 
zweckmäßig zu weiteren Kürzungen zusammengefaßt, nämlich: 


Mit diesen Kürzungen schreiben sich die Ableitungen ö bis yV, die sich von jetzt an auf den Gruppen- 
anfang {=0 beziehen sollen: 


Die Taylorentwicklung der Werte ;=/;A?/2 für i=1 bis 5 vom Ausgangswert 
h/2=}h?j2 aus ergibt: 


11) Siehe z.B. Runge-König:a.a.0.'), S. 287 und Willers:a.a.0O.!), S.305. Der dort verwandte Operator D 
für zwei Veränderliche t, y wird hier sinngemäß erweitert auf drei Veränderliche t, y, v. 
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h? 


+ 12% Dr / +30 DNE +... 
+540 DI +01 DD +... (2). 


2 3 a 
+12 DIE 96 Def + + (640 + 420 Dip 


+ 1440 D}, Di + + 


3 a 
+15 + (150 + 1207, + (1250 + 900 Du 


+ 3000.07, DI +01: +. 3. 


Mit den Werten (11) und (12) lassen sich aus den Gl. (7) und (8) die Fehler 9%, — 9%, bzw. 
h— h anschreiben: 


4 
+ (1205 Di? 2105 Di +5715 DNE +... 


(13). 
+ (140 D? DN 


+ (1125 Di? 1985 Du? 145175 DI + BB +... 


Die Fehler (13) sind jetzt durch eine lineare Verbindung der in Gl. (12) angeschriebenen a-Werte 
nach Ansatz (6) zu Verbesserungen V und V’ derart nachzubilden, daß zwischen Fehlern (13) und Ver- 
besserungen (6) eine Übereinstimmung in möglichst vielen h-Potenzen erzielt wird. Gleichheit in einer 
gewissen A-Potenz verlangt jedoch, wie Gl. (12) und (13) zeigen, von A* an das Abstimmen einer ständig 
wachsenden Anzahl von Einzelgliedern, welche hier bereits als Kürzungen (10) hervorgehoben sind. 
Dem entspricht eine mit fortschreitender A-Potenz wachsende Anzahl von Bedingungen, wodurch die 
Forderung weitergehender Angleichung von korrigierter Näherung und Taylorwert bald ihre Be- 
grenzung findet !?). 


Wünscht man zunächst einen Abgleich bis A*, so sind in 4 und vAh die vier Ausdrücke: 


f für A? 
Di für 
und /„D} für 


zu erfassen. Für die je sechs Beiwerte A; und B; zu den Werten a, bis a, einer Fünfergruppe ergeben 
sich aus der geforderten Übereinstimmung zwischen Fehlern (13) und Verbesserungen (6) in den vier 
Ausdrücken je vier lineare Gleichungen. Je zwei Beiwerte können hiernach frei gewählt, z.B. 
gleich Null gesetzt werden, worauf für die verbleibenden A; und B; die linearen Gleichungssysteme 
heranzuziehen sind. 


6. Ergebnis. Nach Aufschreiben dieser (hier kürzehalber nicht wiedergegebenen) 
linearen Gleichungen zeigt sich jedoch: 


Jedes der beiden Gleichungssysteme ist in sich unverträglich derart, daß der Ab- 
gleich der dritten h-Potenz den gleichzeitigen Abgleich beider Ausdrücke der vierten 
h-Potenz ausschließt. Ein Fehlerabgleich bis h* tritt nur ein, wenn das Glied f, Df, 
d. h. aber die Abhängigkeit der Differentialgleichung von y=v, fehlt. 


12) Entsprechende Überlegungen sind für Differentialgleichungen 1. Ordnung bei der Entwicklung der ver- 
schiedenen Runge-Kutta-Formeln entscheidend gewesen. Vgl. ©. Runge: Math. Ann. Bd. 46 (1895), S. 167 bis 178. — 
K. Heun: Z. Math. Phys. Bd. 45 (1900), S. 23 bis 38. — W. Kutta: Z. Math. Phys. Bd. 46 (1901), S. 435 bis 453. 
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Mit der (Grob-) Korrektur des Blaeß-Verfahrens wird also folgendes erzielt '*): 


1. Für den allgemeinen Fall 5y=f(t,y,j) ergibt die Korrektur volle Überein- 
stimmung zwischen Näherung und Taylorentwicklung bis zur dritten h-Potenz für y 
und vh. Ein Abgleich in höheren Taylorgliedern ist auf diesem Wege grundsätzlich 
nicht mehr möglich. 

2. Für den Sonderfall y=f(t,y), also bei fehlender Abhängigkeit der Differential- 
gleichung von y=v, wird auch das Taylorglied mit y!Y noch voll erfaßt. Ein Ab- 
gleich in weiteren Taylorgliedern ist auch hier nielıt möglich, da die zum nächsten 
Gliede gehörigen Bedingungsgleichungen einander ausschließen. 


Die Zahl der zu einer Gruppe zusammengefaßten Einzelschritte spielt dabei, sobald in 
den a;-Werten alle für den Abgleich wesentlichen Differentialausdrücke erscheinen, keine 
Rolle. Eine Genauigkeitssteigerung kann, wenn erwünscht, mit der in 3. erwähnten Fein- 
korrektur®) herbeigeführt werden. 

Die Genauigkeit der von der (Grob-) Korrektur erzielten Näherung für y und v»h könnte 
hiernach im allgemeinen Falle als „von der Ordnung h°’“ (Fehler von der Ordnung A*), im 
Sonderfalle fehlender Abhängigkeit der Differentialgleichung von % als „von der Ordnung h**“ 
(Fehler von der Ordnung h°) angegeben werden. Zu dieser schon seit Euler üblichen (z. B. 
auch von ©. Runge angewandten) Art der Genauigkeitsangabe ist zu bemerken, daß sie für 
die praktische Genauigkeitsbeurteilung eines Verfahrens nur von bedingtem Wert ist. Bei 
großen Zahlenwerten höherer Ableitungen gewinnen auch die höheren Taylorglieder Einfluß 
auf das Zahlenergebnis, wodurch eine Abschätzung der Fehlergröße auf dem üblichen Wege 
hinfällig wird. — Ein vom Mathematiker angestrebter geschlossener Fehlerausdruck wird sich 
hier, ähnlich wie beim Runge-Kutta-Verfahren, allgemein schwer angeben lassen; auch würde 
er für die praktische Fehlerschätzung nur von geringem Nutzen sein '*). 


7. Eine Korrektur der Zwischenwerte. Näherungswerte im Sinne des Verfahrens von Blaeß sind 
die korrigierten Gruppenend- bzw. -anfangswerte. Schritt im Sinne dieser endgültigen Näherung ist 
also der Gruppenschritt 5A. Die Zwischenwerte 9%, bis 9, und ©, bis , bilden zunächst lediglich 
die Grundlage für den Fehlerabgleich. Sie können indessen in gewissen Grenzen mit Hilfe einer passend 
abgestuften Zwischenkorrektur zur Näherung mit herangezogen werden. Eine solche Korrektur, welche 
die volle Einbeziehung des Taylorgliedes mit % für die Zwischenwerte y; und v; h leistet, ergibt sich 
auf Grund der Angaben in 5. Sie lautet: 


1 
vh=uh+—V 
5 


12 (da). 


II. Verfahren nach dem Runge-Kutta-Prinzip. 


8. Runge-Kutta-Näherung für Differentialgleichungen 2. Ordnung. Es erscheint lohnend, 
die im Blaeß-Verfahren für Differentialgleichungen 2. Ordnung angewandte Rechnungsform 
auf eine Methode zu übertragen, welche bei Differentialgleichungen 1. Ordnung in ähnlicher 
Problemstellung zum Erfolg führt. Wir meinen das Prinzip der Runge-Kutta-Näherung. 

Das Runge-Kutta-Verfahren zur Lösung von Differentialgleichungen 1. Ordnung") be- 
dient sich verschiedener vorläufiger Steigungswerte ij, oder besser vorläufiger Tangenten- 
erhebungen k;=yh, genommen an verschiedenen Stellen im Schrittinnern. Aus diesen 
Werten k; wird ein Mittelwert k so gebildet, daß der damit berechnete Neuwert y,=9y,+ k 
am Schrittende in möglichst vielen Gliedern mit einer Taylorentwicklung übereinstimmt. 

Überträgt man dieses Prinzip der Mittelung auf Differentialgleichungen 2. Ordnung, so 
treten an die Stelle der Steigungen 5 die Wölbungen j, an die Stelle der Erhebungen 


13) Für den einfachen Sonderfall y —=/(t) erreicht die Korrektur nach GI. (4) nieht nur vollen Fehlerabgleich 
in 4% für y und oh, sondern, herbeigeführt durch die getroffene Auswahl der drei «-Werte für V und F’ aus den 
Werten «@o bis a,, eine weitgehende Annäherung in den höheren Taylorgliedern weit über y!Y hinaus. Für diesen 
Sonderfall y=/(t), also für den Fall der zweifachen mechanischen Quadratvr, erzielt schon die (Grob-)Korrektur so 
hohe Genauigkeiten, wie sie im allgemeinen Falle die Feinkorrektur herbeiführt. 

14) Vgl. hierzu R. v. Mises: Zur mechanischen Quadratur. Z. angew. Math. Mech. Bd. 13 (1933), S. 53 bis 56. 

15) Darstellungen siehe a. a. O.t). 
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k;=Yh die Abwölbungen a;=ih?]2. Aus vorläufigen 

a;-Werten, berechnet an Zwischenstellen {,+4h im Schritt- 

innern mit Zwischenwerten y und v, sind zwei Mittel- 


x werte a und a’ zu bilden, aus denen die beiden Neuwerte 
4 y, und v,h in ,=t,+h berechnet werden können in einer 
„| der Gl. (8) ähnlichen Form: 

y-ytvh+ta, .. (14). 
| v Diese Endrechnung wird in Bild 3 veranschaulicht. 

j 


& 


Die Zwischenwerte y und v»h werden nach Schritt- 
. gleichungen errechnet, die aus (3) durch Umschreiben auf den 
1 % Teilschritt Ah entstehen, wobei grundsätzlich alle bis dahin 
im Schritt bekannten Werte a; in geeigneten Anteilen o;a; 
to "2 bzw. 7;a; herangezogen werden können. Die Zwischenwerte 


. sind also von der Form: 
Bild 3. Aufbau der Näherungswerte h +2: 


nach GI. (14). 
ceh=v,h+2i 


Für die Zwischenstellen werden zweckmäßig die in der gebräuchlichen Runge-Kutta-Formel 
bewährten Bruchteile 


1 
vorgesehen. Aus den an diesen Stellen berechneten vier Werten a,, a,, a,, a, werden die 
zwei Mittelwerte « und a’ zur Berechnung der endgültigen Näherungen 4, und v,h am 
Schrittende aufgebaut. 

Das Ergebnis einer nach diesem Plan durchgeführten Rechnung ist nicht neu. Die ge- 
wonnenen Formeln stimmen inhaltlich mit Formelsätzen von Nyström?°) überein, die dieser 
— in Forın und Herleitung anders — für Differentialgleichungen 2. Ordnung und für Systeme 
solcher Gleichungen angegeben hat. Der hier verfolgte Weg zeigt die Formeln nach Auf- 
fassung und Handhabung in neuem Licht: 

1. Die Nyström-Formeln erscheinen hier als folgerichtige Fortentwicklung der 
Runge-Kutta-Näherung auf Differentialgleichungen 2. Ordnung, in- 
dem die der Mittelung unterzogene Wölbung % an die Stelle der die Differential- 
gleichungen 1. Ordnung kennzeichnenden Steigung j tritt. Deı Gedanke Nyströms, 
das zuvor übliche Vorgehen des Aufspaltens in zwei Differentialgleichungen 1. Ord- 
nung'*) durch unmittelbare Behandlung der Differentialgleichung 2. Ordnung zu er- 
setzen '”), wird durch unsere Betrachtungsweise anschaulich bestätigt. 

2. In der hier gewählten Form erhalten die Nyström-Formeln eine Fassung, die das 
Verfahren namentlich auch dem wissenschaftlich arbeitenden Ingenieur als ein zu- 
verlässiges und in der Handhabung einfaches Instrument zur Lösung von Differential- 
gleichungen 2. Ordnung zugänglich macht. 


9. Neue Fassung der Nyström-Formeln. Die Durchführung eines Nyström-Schrittes h, 
der einer Schrittgruppe 5h des Blaeß-Verfahrens ähnlich ist, gliedert sich wie folgt: 

1. Vorbereitende Schrittrechnung zur Gewinnung von vier vorläufigen Werten 
a; an den Zwischenstellen t— t,—=0, h/2, h/2, h: 


h? 

a=fl(b; 
h 1 1 h? 


w 


h 1 1 


h? 
a=fl,+h; ytuht+a,; 


An Stelle der Werte y=v sind die mit y dimensionsgleichen Werte vh verwendet. 
2. Mittelbildungen a und a’: 


16) Siehe z. B. Fr. A. Willers, a. a. O.1), S. 335, Gl. (65). 
17) Nyström,a a. 0.5), 8.7,$2. 
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3. Endreehnung der Neuwerte: 
y=y+tv,h+ta, 


Die Rechnung wird übersichtlich in folgendem Schema nach Tafel 2 durchgeführt (beginnend 
mit £,—=0): 


Tafel 2. Nyström-Rechnung in neuer Fassung. 


t vh | @; a,a’ 
0 h a; a 
| 1 1 
| 4 a, | vh+a, a’ 
a, 
ah 
| 


Zunächst wird Zeile für Zeile eines Schrittes bis zum Doppelstrich gerechnet. Es 
wechseln ab: Bestimmung des a;-Wertes durch Einsetzen der in gleicher Zeile stehenden 
Werte i,9,vh in die Differentialgleichung; Berechnung der Zwischenwerte y und vh der 
neuen Zeile mit dem voranbestimmten a;-Wert in der im Schema angegebenen Weise. — 
Nach Vorliegen des letzten a-Wertes a, erfolgt hinter dem Doppelstrich die Mittelbildung «a 
und a’ nach (16); danach die Endrechnung der Neuwerte nach (14). 

Die Werte y; und v;h der ersten Zeile jedes Schrittes werden als die endgültigen 
Näherungen zweckmäßig unterstrichen, da auf sie bei der Aufstellung der Zwischenwerte % 
und vh der drei folgenden Zeilen und bei der Endrechnung zurückgegriffen wird. Die Be- 
rechnung des y-Wertes der beiden Mittelzeilen ist bei einer ohnehin zu empfehlenden Ma- 
schinenrechnung bequem durch Eintasten der ganzen Werte v,h und a, und Multiplizieren 
mit 0,5 bzw. 0,25 auszuführen. 


Für Sonderfälle, in denen % unabhängig ist von y oder , ergeben sich willkommene 
Vereinfachungen °®): 


A. Sonderfall: y=f(t,v), unabhängig von y. 

Die Zwischenwerte für 4 brauchen nicht berechnet und im Schema mitgeführt zu 
werden, wodurch sich die Rechnung vereinfacht. Mittelbildungen (16) und Endrechnung (14) 
bleiben unverändert. Inhaltlich stimmt die Rechnung jetzt mit den Runge-Kutta-Formeln für 
zwei Differentialgleichungen 1. Ordnung überein. Erläuterung durch Zahlenbeispiel 1, Tafel 4. 


B. Sonderfall: y=f(t,y), unabhängig von y=v"?). 

Die Zwischenwerte für v»h brauchen nicht berechnet und im Schema mitgeführt zu 
werden. Da die - und y-Werte der beiden Mittelzeilen je in sich gleich sind, so verkürzt 
sich das Schema um eine Zeile auf nur drei Zeilen, ein gegenüber den alten Runge- 
Kutta-Formeln'*) bei der Wichtigkeit dieses Sonderfalles wesentlicher Vorteil. Mit ab- 
geänderten Bezeichnungen (a,, a,, a,) lautet das Schema: 


Tafel 3. Nyström-Rechnung in neuer Fassung für y=f (t, Y). 


t | Yy a; a,a 
1 1 ‚ 
hj2 | @ 
h as | 
h Yı | 
| 


Für die Mittelbildung ist der Wert a, der mittleren Zeile doppelt zu zählen: 


1 
a=7(a +2a,), a 
Die Endrechnung erfolgt unverändert nach Gl. (14). 


18) Der Sonderfall y=/(t) führt auf eine der Simpsonregel entsprechende Formel für zweifache Quadratur, 
19) Dieser Fall entspricht den Nyström-Formeln Nr. II, Nyström, a. a. 0.5), S. 24 
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Für Systeme von Differentialgleichungen 2. Ordnung werden die an- 
gegebenen Formeln nach Nyström entsprechend erweitert. Das Rechnungsschema ergibt 
sich durch Vervielfachung der in Tafel 2 bzw. Tafel 3 dargestellten Anordnung gemäß 
mehreren Veränderlichen &, %,... . Namentlich die in der Technik wichtigen Differential- 
gleiehungen 4. Ordnung lassen sich so durch Aufspalten in ein System von zwei 
Differentialgleichungen 2. Ordnung behandeln. Eine unmittelbare Behandlung erscheint nicht 
zweckmäßig, weil im Rechenschema unbequeme Vorzahlen auftreten würden. 

Die Genauigkeit des Nyström-Verfahrens®) ist ähnlich der des Runge-Kutta- 
Verfahrens von der Ordnung h* (Fehler von der Ordnung h°) für y und v. Ein praktisch 
brauchbarer Fehlerüberschlag, der sich auch zur Verbesserung der Werte y (und 
ebenso von v» bzw. vh) benutzen läßt, ergibt sich aus einer mit doppelter Schrittweite 2h 
wiederholten Rechnung (Grobrechnung). Die allerdings nur bei genügend kleiner Schritt- 
größe h sinnvolle Regel lautet?'): 


1 
d= — =y+ 4 (17). 


Die Wahl der Schrittgröße Ah richtet sich nach der geforderten Genauigkeit; eine 
allgemeine Regel läßt sich dafür nicht angeben. Die Genauigkeit ist i. a. um so größer, je 
weniger der letzte Wert a, (a, im Sonderfall B) des alten Schrittes vom «,-Wert des neuen 
Schrittes verschieden ist. Der Schritt ist sicher zu groß, wenn beide a-Werte stark von- 
einander abweichen. Die Schrittgrößen h sind u. U. im Laufe der Rechnung entsprechend 
abzustufen. 

Zur Erläuterung des Rechnungsganges folgen zwei einfache Zahlenbeispiele, und zwar 
für die Sonderfälle A und B. Die Durchrechnung eines allgemeinen Falles wird in 11. 
wiedergegeben. 


1. Beispiel: Strenge Lösung: y=1-— e-?t 


Y%= t=0,1: 0,181 269 

vh— 0,163 746 
Schritt: h=0,1 


Tafel 4. Beispiel y= — 2%. 


02 0980067 | — 0,039 734 


t | Y | vh a; a,a’ 
0,00 0 | 0,200 000 — 0,020. 000 — 0,018 733 
0 | 180000 — 18000, — 17520 
0 | 182000 | — 18200 
10 | 163600 °— 16360 | 
0,10 | 0,181267 | 0,163747 | 
| 
2. Beispiel: — Strenge Lösung: y=cost 
Y = t=0,2: Y= 0,980 067 
vh= — 0,039 734 
Schritt: h=02 
0,02 
Tafel 5. Beispiel — y. 
t Yy oh | a; | a,a 
0,0 1,000 000 | 0 F% 0,020 000 | — 0,019 933 
1 995000 | | — 19801 
980 100 | I— 1902| 
| 
| 


20) Nyström, a. a. 0.5), S. 26 bis 27. Vgl. auch das hierzu grundsätzlich unter 6. Gesagte. 
21) Vgl. z. B. Willers, a. a. O.!), S. 310 oder auch Nyström, a. a. 0.5), S. 27. 
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10. Näherungsverfahren für Differentialgleichungen 3. Ordnung. Das Runge-Kutta-Prinzip 
läßt sich sinngemäß auf Differentialgleichungen 3. Ordnung erweitern. Der Schrittrechnung 
liegt dann die nach ij abgebrochene Taylorformel 


zugrunde. Die Glieder 4,, =r, und 7, sind als Anfangswerte gegeben, das Glied mit 
wird durch Einsetzen der Anfangswerte in die Differentialgleichung 


y=fl,yv,w) 


berechnet. Für das Glied ö h?/6 ist eine Mittelbildung in ähnlicher Weise vorzunehmen, wie 
für Jh beim Runge-Kutta-Verfahren für Differentialgleichungen 1. Ordnung oder für die 
a;-Werte beim Nyström-Verfahren. 
Mit Rücksicht auf eine bequeme Zahlenrechnung wird der vom Taylorglied abweichende 
Ausdruck 


benutzt. Behält man die bewährte Schritteinteilung t=0,h/2, h/2,h bei, läßt jedoch beim 
Aufbau der Zwischenwerte für y, v und w zunächst noch gewisse Verallgemeinerungen 
offen”), so führt die Durchrechnung unter Voraussetzung eines Abgleichs bis zur Ordnung 
h* in 9, v und w auf einen Formelsatz, welcher einer Aufspaltung der Gleichung 3. Ordnung 
in drei Gleichungen 1. Ordnung und Durchrechnung nach den alten Runge-Kutta-Formeln 
hinsichtlich Arbeitsaufwand und Genauigkeit überlegen ist. Für die beiden Mittelzeilen t=hj2 
sind, ähnlich wie bei den Nyström-Formeln, diesmal die Werte y und oh je in sich gleich. 
Hieraus und aus der einer Taylorformel angepaßten Anordnung ergibt sich eine namentlich 
in Sonderfällen fühlbare Arbeitsersparnis, welche den zunächst etwas verwickelteren Aufbau 
der Schrittrechnung wettmacht. 


Der Rechnungsgang verläuft nach dem Schema Tafel 6 ähnlich dem für die Nyström- 
Formeln angegebenen, aus welchem die Schrittrechnung genügend klar hervorgehen dürfte, 
und nach den Gl. (20) und (21) für Mittelbildungen und Endrechnung. Die Berechnung der 
Mittelzeilen für y läßt sich bei einer ohnehin notwendigen Maschinenrechnung bequem durch 
Eintasten der Werte v,h, a, und £,/3 und Multiplizieren mit den Faktoren 0,5, 0,25 und 0,25 
ausführen; in ähnlicher Weise in der letzten Zeile £,/3 mal 2. 


Schrittreehnung: 


Tafel 6. Rechnungsgang für eine Differentialgleichung 3. Ordnung. 


t | Yy ch a 
uch 
| | 
h/2 | ” ” ” ” + | 


22) Der Rechnung der beiden Mittelzeilen des Schemas Tafel 6 liegt folgender allgemeiner Ansatz zugrunde: 
1 1 1 
1 
vh=voh+%+ zpı 
+ Pı 
Zeile III: y=yo+ voh +4 0+ 5 +uPße] 
1 
3le-9 Bı+ 


Pr]. 


Nach der vereinfachenden Annahme $=1 ergeben sich die unbestimmten Größen 4, «, ©, r, » in eindeutiger Weise aus 
den gestellten Genauigkeitsforderungen für 4, v und w. 


8 
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1 1 1 
Mittelbildung®): 


Endreehnung: 


Für die Sonderfälle, in denen die Differentialgleichung unabhängig ist von y oder %, 
brauchen die Zwischenwerte für y bzw. vh nicht berechnet und im Schema mitgeführt zu 
werden (1. Beispiel). 

Für den Sonderfall fehlender Abhängigkeit von y verkürzt sich das 
Schema ähnlich wie beim Nyström-Verfahren auf nur drei Zeilen. Der Vorteil gegen- 
über der alten Runge-Kutta-Rechnung ist hier besonders augenfällig. Die Rechnung verläuft 
wie folgt: 

Schrittreehnung: 

Tafel7. Rechnungsgang für Y=f(t, y,Y). 


t Y ch a 


0 Yo h 


vh+20+28, 


h a, 
| | 
Mittelbildung: +B) 


Endrechnung wie oben nach Gl. (21). 
Zur Erläuterung folgen zwei einfache Beispiele: 


1. Beispiel: Yy— — Strenge Lösung: 
Y=%=V0 y=taretgt— In(1+) 
1 _ 
— a = 0,019 2308 
Tafel8. Beispiel 
0,0 0.020 0000 0 — 1307 
1 | | 20 0000 — 0,000 4000 — 2614 
1 | 19 6000 3842 — 5080 
DAR | 19 2317 7397 
02 | 0,0198693 | 0,039 4772 | 0,019 2307 
| 


23) Mit einer Mittelbildung ?=1B0-(9 +6ß2+6ßs— wird für y sogar ein vollständiger Abgleich des 
Gliedes mit A5 erzielt. Hier ist die zahlenmäßig bequemere Mittelbildung nach (20) vorgezogen, die sich der genauen 
Formel gut anpaßt. 
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2. Beispiel: Y—=y Strenge Lösung: y=el 
t= 0,2 1,221 4028 
v h— 0,244 3806 
Schritt: h=02 a — 0.024 4281 
8; — 0.002 y 
Tafel 9. Beispiel 
y oh a | 
0,0. 1,000.0000 | 0,200. 0000 | 0,020.0000 0.0020000 14034 
1. 1,105 1667 | | 22103 2 1402 
2 | 1,221 4735 24429 2 2878 


02 | 12214034 | 0,244 2804 


| 


0.024 4280 


| 
| 


Zum Vergleich der Näherung mit den genauen Taylorwerten für y, vh und a zeigt 
Tafel 10 die Faktoren der einzelnen Differentialausdrücke in den Gliedern mit h° und h®, be- 
zogen auf den genauen Wert 1, und zwar in Gegenüberstellung zum alten Runge-Kutta-Ver- 
fahren mit Aufspaltung in drei Gleichungen 1. Ordnung. Der Operator D ist, ähnlich wie 
unter 6, als Kürzung für eine vollständige Differentiation nach ? anzusehen, jedoch hier für 
Funktionen von t#, y, v, w. . 


Tafel 10. Vergleich der Näherung mit genauen Taylorwerten 1 
in den Potenzen A® und A®. 


Y ch a 
h- Faktor 
Dobias bei Neue Runge- | Neue | Runge- | Neue , Runge- 
Formel | Kutta | Formel | Kutta | Formel | Kutta 
54 5/4 1 1 1 1 
h5 il» 5/4 0 1 | 1 1 1 
h5 Dite 5/4 0 1 1 1 1 
hs Ds} 5/4 5/4 5/6 5/6 1 1 
hs 3w 5/4 0 5/6 5/6 
hs 3/D}, 4 | 0 5/6 5/6 1 1 
hs 5/6 5/6 1 1 
hs fly 5/4 | 0 5/6 0 1 1 
hs Di}, 0 0 0 0 1 1 
hs 15/8 0 5/4 5/4 1 1 
hs wfyhe 15/8 0 5/4 0 1 1 
h® Hals 15/8 | 0 5/4 0 1 1 
hs 0 | 0 1 1 


11. Durchführung eines Beispiels. Das Verfahren von Nyström wurde zur numerischen Integration 
einer der Besselschen Differentialgleichung verwandten inhomogenen Gleichung angewandt, nämlich: 


1 
mit den Anfangswerten = — 1 für wofür R. Müller?*) Zahlenwerte mit Hilfe mecha- 


nischer Quadraturen berechnet hat. 


Die numerische Integration erfolgt dabei zwischen Kontrollwerten (in Tafel 12 mit * versehen), 
welche durch Reihenentwicklung berechnet worden sind. Das Verfahren arbeitet also hier nach Art 
einer Interpolation. Die Genauigkeit wurde in der unter 9. beschriebenen Weise nach Gl. (17) gesteigert. 
Die Fehler betragen bis zu zwei Einheiten der letzten in Tafel 12 angegebenen Stelle. 


24) Reinhard Müller: Zahlenmäßige Beherrschung und Anwendung Besselscher und verwandter Funktionen. 
Z. angew. Math. Mech. Bd. 19 (1939), S. 36 bis 54. — Das Beispiel wurde auf Anregung von Herrn Professor 
Dr. ©. Schmieden, Darmstadt, durchgeführt. 
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Tafel 11. Anordnung der Rechnung. 


t Yy | ch a; | a, a 
1.00 0,8730 8424 | — 0,0468 4508 | 0,0017 0768 16 6229 
05 8500 8862 451 3740 16 3794 16 1775 
05 s | 452 0714 16 4126 
10 8278 8042 | 435 6256 15 7406 
1,10 0,8279 0145 | — 0,04356504  0,00157428 | 153268 
15 8065 1250 | 419 9076 15 1042 14 9177 
| 420 5462 15 1320 
20 TBB 405 3864 14 5169 
1,20 0,7858 6904 | — 0,0405 4064 
Tafel 12. Ergebnisse für das Beispiel d+ 
t Y t | Y | re 
0,0 1,5707 9633 1.0000 0000 40. *026810472 |  *0,0729 353 
1 * 4746 1462 * 0,9246 9752 de 2545193 660 678 
2 * 3856 5329 * 8555 0259 24921 &00 326 
3 3033 279 7919 003 6. 2304584 547 152 
4 * 22710253 * 7334 2035 8 | 2199 947 500 178 
5 1564 873 6796 329 50 | * 21041555 | * 4585690 
6 0910 327 6301 447 8 | 2016206 421610 
7 1.0303 270 5845 970 4 | 1935245 388 697 
8 0,9739 929 5426 613 6 | 18604963 359 303 
9 9216 844 5040 374 8 1791 3158 332 981 
1.0 * 0,8730 8424 * 0,4684 508 6.0 0,1727 1246 0,0309 349 
1 8279 014 4356 504 2 16674198 288 073 
2 7858 690 4054 064 4 11155 268 868 
3 746749 | 3775 087 6 1559 7504 251489 
4 7102955 | 3517646 8 1511 0546 235 721 
5 0763231 | 3279 984 7.0 1465 3670 221379 
6 6446352 3060 488 2 1422 4185 208 303 
7 615057 2857 682 4 1381 9710 196 351 
8 | 58743065 | 2670 214 6 1343 8113 185 403 
| 2496 847 8 1307 7500 175 351 
2,0 * 0,5374 5039 * 0,2336 445 8,0 0,1273617° 0,0166 102 
2 4936 545 2050 468 2 24161 157 574 
4 4551 621 1804 965 #: 1210 548 149 693 
6 4212 272 1593 782 6 1181 345 142 398 
8 3912 167 1411 740 8 1153 551 135 631 
3,0 3645 926 1254 477 9.0 1127 061 129 342 
2 3408 971 1118315 2 1101 784 123 489 
4 3197 401 1000 149 4 1077 639 118 036 
6 3007 887 897 358 6 1054 548 112 932 
8 2837 582 807 721 8 1032 444 108 162 
we 10,0 * 0,1011 2644 * 0,0103 6927 


12. Zusammenfassung. Nach einer kurzen Darstellung des Verfahrens von Blaeß 
zur angenäherten Lösung von Differentialgleichungen 2. Ordnung werden einige die (Grob-) 
Korrektur des Verfahrens betreffende Fragen auf dem Wege von Taylorentwicklungen unter- 
sucht. Über den erreichbaren Grad der Übereinstimmung mit einer Taylorreihe lassen sich 
bestimmte Aussagen machen. 

Zu der von Nyström herrührenden Übertragung des Runge-Kutta-Verfahrens 
auf Differentialgleichungen 2. Ordnung wird, aufbauend auf der Blaeßschen Rechnungsform, 
eine handliche Fassung angegeben, der sich auch Sonderfälle leicht einfügen. 

Auch für Differentialgleichungen 3. Ordnung lassen sich auf der Grundlage 
des Runge-Kutta-Prinzips Formelsätze herleiten, welche eine unmittelbare Behandlung dieser 
Gleichungen gestatten. Gegenüber dem Aufspalten in drei Differentialgleichungen 1. Ordnung 
ergeben sich, namentlich in Sonderfällen, gewisse Vereinfachungen. 

Die Anwendung der Rechenverfahren wird durch einige Beispiele erläutert. 132 
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Zur Umwandlung von rotationssymme- 
trischen Problemen in zweidimensionale 
und umgekehrt. Auf eine Halbebene 4<0 
wirke die rotationssymmetrische Belastung p, (r). 

Beispiel: (Bild 1) ()=cfürr<a, p, (r)=0 
für r >a. 

Wir bilden hieraus einen zweidimensionalen Be- 
lastungsfall durch einen Translationsprozeß T, in- 
dem die Belastung in z-Richtung verschoben und 
auf die Länge 2 verteilt wird. Die neue Belastung 
wird: 


+» 
f 


Für das Beispiel wird 9%, =cyYa? — 2? für | <a 
und für |z|)>a (Bild 2). 

Nun vollführen wir die Operation D, die in 
folgendem besteht: Wir nehmen p, als Randwerte 
einer analytischen Funktion der Halbebene y<0, 
die für |e+iy|>» verschwindet, und differen- 
tieren sie am Rande 0) nach Die so erhaltene 
Funktion ist: 


Ps (X) = 


R f». 


| 


0 


c Bild 1. 
l (links oben) 


a 


Bild 2. 
(links unten) 


Bild 3. 
(rechts unten) 


1 


Für unser Beispiel wird die analytische Funktion 


(tip? —ilctiy). 


Die Ableitung nach wird eR (7: + 


Am oberen Rande der Halbebene erhält man (Bild 3) 


ilz 


Ps(2)=c für und Vz) 
für |z|>a. 

Durch einen Rotationsprozeß R bilden wir aus 
dem zweidimensionalen Belastungsfall p, (z) einen 
rotationssymmetrischen 2, (r). Der Prozeß besteht 


darin, daß man den Belastungsfall 9, um die 4-Achse 
um den Winkel dreht und jeweils den Teil 


ne wirken läßt. Im Grunde genommen bilden 


wir für jede Kreislinie mit dem Halbmesser r um 
die y-Achse den Mittelwert von 2, (2). 


zn 


Für unser Beispiel erhält man 
für r<a und für r>a. 


Wir haben die Ausgangsfunktion p, (r) zurück- 
erhalten. Es folgt für das Beispiel: 


Dieses gilt aber auch für jede andere rota- 
tionssymmetrische Belastung p, (r), die als Über- 
lagerung endlich- oder unendlichvieler als Beispiel 
behandelter Belastungsfälle gebildet werden kann. 

Diese Prozesse lassen sich mit allen Größen vor- 
nehmen, die sich wie Skalare überlagern, also für 
die Verschiebungen ® in Richtung der Rotations- 
achse bei elastischen Problemen, für oy, &y, 
s=0,+0y+0%, +8; bzw. 


rdr 


1 
Rp, (2) = 


ur. 


| Y 


— a 
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Als Anwendungsbeispiel werde gezeigt. wie die 
Spannungen an der Oberfläche und auf der Achse 
beim Belastungsfall nach Bild 1 gefunden werden 
können, da der Spannungszustand nach Bild 2 
bekannt ist. 

Für den zweidimensionalen Spannungszustand 
der Belastung 9, (7) wird, mit 


=2ug,, 
1 
s=2(1+ug9, &= 


Für den zweidimensionalen Spannungszustand 
der Belastung p (2) nach Bild 3 ist statt 9, zu 


setzen: 
i(ce-+iy 
= cRi1i— 
| 
An der Oberfläche „=0 wird für p;: 
_ 
E Ps, 
2(1 1—2u) 


E 


Für den rotationssymmetrischen Fall nach Bild 1 
erhält man durch den Prozeß R: 


oy=—c für r<a bzw 0 für r>a, 
bzw. 0, 


= — 
E bzw. 0. 
Hieraus für r<a 
bzw. für r>a, 
_du 
2E 
a 
„4 
r 2E 
a: 
bzw. — SE 5 
a: 
3 bzw. + — 3 


2 2 r?" 
Auf der Achse wird für p, sowohl wie für p, 


hierbei ist „< 0 


Damit ist auch o,—=0; auf der Achse des rotations- 
symmetrischen Falles nach Bild 1 gefunden. 


Dresden. C. Weber. 161 


Die erste Variation der Formänderungs- 
arbeit ausgebeulter ebener Platten. 


Bei Behandlung von Problemen der Platten- 
bierung, Plattenschwingune und Plattenkniekung 
geht man am besten davon aus, daß die potentielle 
Energie bzw. das Hamiltonsche Integral einen Ex- 
tremwert haben muß. Denn aus der Variation 
dieser Ausdrücke folgt nicht nur die Differential- 
gleichung, sondern auch sämtliche Randbedingun- 
gen und bei sprungweiser Änderung der Platten- 
dicke oder bei Einlagerung linear verteilter, mit- 
schwingender Massen auch die Übergangsbedingun- 
gen. Sowohl in den Ausdruck für die potentielle 
Energie (Differenz der Formänderungsarbeit und der 
Arbeit der äußeren Kräfte) wie in den Integranden 
des Hamiltonschen Integrals (Differenz der poten- 
tiellen und der kinetischen Energie) geht die Form- 
änderungsarbeit ein. Der Ausdruck für diese ist 
z. B. von Kirchhoff mittels der Hauptkrünmungen 


und der ausgebeulten Mittelfläche der Platte, 


d. h. der Fläche, die vor der Deformation in der 
z, y-Ebene lag, gegeben: 


+ 


wo N der Plattenfaktor, und nl die Querkon- 


traktionszahl ist. Für kleine Dinelliehnäeen geht 
das, wenn man mit w (z, %) Verschiebungen senk- 
recht zur 7, Y-Ebene bezeichnet, in 


über. Vorausgesetzt ist, daß jede Gerade, die ur- 
sprünglich senkrecht zur Plattenmittelebene stand, 
durch die Deformation in eine Gerade senkrecht zur 
verbogenen Plattenmittelebene übergeht, und daßalle 
Teile der Mittelebene beim Ausbeulen keine Deh- 
nung erfahren. Die Doppelintegrale sind dabei über 
das ganze Plattengebiet zu erstrecken. Während 
nun die Variation des ersten Integrales von A sich 
einfach mittels des Greenschen Satzes behandeln 
läßt, fordert die Behandlung des zweiten 


das keinen Einfluß auf die Differentialgleichung hat, 
sondern nur auf Übergangs- und Randbedingungen, 
und das, wie man leicht feststellt, der Form nach 


invariant gegen Drehungen des Koordinatensystems 
ist, eine etwas umständlichere Umformung. 


1. Die erste Variation von J ergibt: 

dy? 
"wd?öw 

- Jazas 


droydrdy 
Hieraus erhält man durch unvollständige Integration 
des ersten Summanden im ersten Integral nach x, 
des zweiten nach Y 


(1). 


r o (2), 
u u 
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da das Flächenintegral auf der rechten Seite Null 
wird. Beachtet man nun, daß nach Bild 1 


dy=dssind, de—=dscos» 
ist, so kann man dieses Integral in eine Summe von 


Integralen über die einzelnen Bogen des Randes von 
Ecke zu Ecke zerlegen, wozu man, wenn man den 


€ 
& 
X 
&, 
& 
dy 
7 
€ 


Bild 1. 


Umlaufssinn im Gegenzeigersinn wählt, das Vor- 
zeichen des zweiten Integrales umkehren muß. Es 


wird 
Ad sind 
(3), 
+ 2. — .. 
drdy 


wo &, die Bogenlänge bis zur »-ten Ecke gemessen 
von einem beliebig gewählten Anfangspunkt aus 
ist; dabei sind die Koordinaten von *, gleich denen 
von &,. In der gleichen Weise formt man das 
zweite Integral von (1) um. 


2 
‚=1 (4). 


als 
Zur weiteren Umformung beachtet man, daß, falls 
n die äußere Normale ist, 


Dow. 
35 cosd—+ sin 


Setzt man das in die Integrale (3) und (4) ein, so 
wird 


s=l &-ı+ 


sind + — —sin?®d) 
6). 
”w ., 


w 
+2, din ds 
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Da nun 
dy®dn 


ist, wird der Ausdruck in der zweiten geschweiften 
Klammer 


dw 


Ähnlich findet man durch Differentiation nach s, 
wobei man berücksichtigen muß, daß #, das ja von n 
unabhängig war, sich mit s ändert: 


Im dsind (cos? — sin? ®) 
w dd 


Da es sich um Kurvenbogen ohne Ecken handelt, 
kann man für z die Krümmung des Randes > 
einsetzen, ferner ist 


dw dw 
sin 


so daß also die erste Klammer des Integrales (5) wird 


Setzt man (6) und (7) in (5) ein, erhält man 


dw 


— N 


dw 
w 1dw 
+ 


Das erste dieser beiden Integrale wird nochmal 
durch unvollständige Integration umgeformt. Dabei 
sei angenommen, daß in m Punkten o, des Randes 
Einzelkräfte P,, senkrecht zur Plattte wirken. Der 
einfacheren Schreibweise wegen sollen jetzt auch 
die Ecken mit o, bezeichnet werden und die m+n 
=r Punkte, in denen Einzelkräfte angreifen bzw. 
in denen der Rand eine Ecke hat (natürlich können 
auch in diesen Ecken Einzelkräfte angreifen), fort- 
laufend numeriert werden. Dann wird, wenn man 
im ersten Gliede auf der rechten Seite die Reihen- 
folge der Summanden etwas umordnet: 


dw 
dsdn ds 
dw, 1dw"+" 


Setzt man jetzt wieder statt der Summen von Inte- 
nt Integrale über den ganzen Rand, indem man 
nstetigkeiten des Integranden an einer endlichen 
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Anzahl von Punkten zuläßt, erhält man als End- 
resultat 


r 
® dw 1 d +0 


dw »/1dw . (8). 


+ 

2. Plattenbiegung. 

Die potentielle Energie, die ja gleich der Differenz 
aus Formänderungsarbeit und Arbeit der äußeren 
Kräfte ist, läßt sich für die Plattenbiegung, wenn 
außer den oben erwähnten Einzelkräften am Rande 
noch eine Flächenbelastung » (x, y) und eine konti- 
nuierlich verteilte Randkraft P (s), beide positiv in 
Richtung positiver w und außerdem ein Rand- 


dw 
moment M (s), das das gleiche Vorzeichen wie = 


hat, an der Platte angreifen, schreiben: 


u=5 [[uwraray 


- 
vn 
r 
- Pro Po, 


Bildet man die erste Variation dieses Ausdruckes, 
wobei angenommen sei, daß der Rand so beschaffen 
ist, daß man für die Umformung der Variation des 
ersten Integrales die Greensche Formel benutzen 
kann, so wird 


öwddwdrdy 


+ 
r 


osdn 
dw d/1dw 
+ 9). 
1dw 


4 
- wo, = 0 


Soll diese Gleichung erfüllt sein, muß, da dw 
innerhalb der Begrenzung der Platte beliebig sein 
kann, 


NAAw— p(e,y)=V0 


sein. Damit hat man die Differentialgleichung der 
Plattenbiegung, in die nichts von dem Integral J 
eingeht. Das Nullsetzen der anderen Glieder gibt 


die Randbedingungen. 
w 


on 
muß, wenn (9) erfüllt sein soll, 


Bei frei drehbarem Rand, 


ganz beliebig genommen werden kann, 


wo also 


Hier rührt der zweite Summand von J her. Ist 
dagegen der Rand eingespannt, so werden die ent- 


sprechenden Glieder Null, weil ne gleich Null zu 


setzen ist. An die Stelle der Randbedingung Ia 
tritt 
lb: 


und aus Ia berechnet sich dann, wie man sich leicht 
überlegt, das Auflagermoment. 

Bei frei verschiebbarem Rand, wo dw(s) und 
öw, beliebig sein können, muß zur Erfüllung der 
Gleichung (9) gesetzt werden 


E ” dw /1dw 


Hier kommt in der ersten Gleichung der zweite 
Summand rechts von J her, während die zweite 
Gleichung überhaupt nur von J herrührt. Ist die 
Platte drehbar gelagert, so wird öw(s)—0 und 
öWs — 0, an Stelle von IIa tritt 


Ib w(s)=0 und w=0 


und aus Ila berechnen sich die Auflagerreaktionen. 
Natürlich werden jetzt im allgemeinen keine Einzel- 
kräfte auftreten, da in diesem Falle der Wert der 
eckigen Klammer für o +0 und o—0 der gleiche 
sein wird. Eine Ausnahme machen nur die Ecken. 
Hier werden Einzelkräfte auftreten: 


oder, wenn das nicht der Fall sein kann, ist durch 
das Nullsetzen der rechten Seite eine Bedingung für 
w in dieser Ecke gegeben. Diese sog. Eckenreiation 
rührt also nur von J her. 

Im Fall der Plattenschwingung, auf den hier nicht 
weiter eingegangen werden soll, hat man, falls keine 
mitschwingenden Einzelmassen vorhanden sind, die 
gleichen Randbedingungen, wie bei der Platten- 
biegung. Nur in der Differentialgleichung tritt noch 
ein weiteres Glied mit der zweiten Ableitung von w 
nach der Zeit auf. 


3. Plattenknickung. 


Die potentielle Energie ist jetzt durch die sog. 
Bryansche Formel gegeben: 


y? 
w\? 

wo P, und P, die in der Plattenmittelebene auf die 
Längeneinheit wirkenden Normaldrucke in der x- 
bzw. der y-Richtung sind, und wo T die Schub- 
spannung pro Längeneinheit ist. Die Variation 
dieses Integrales soll hier nicht hingeschrieben 


2 
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werden. Es soll nur kurz der Einfluß von J auf die 
Beullast betrachtet werden. Diese bestimmt sich 
bekanntlich durch Nullsetzen einer Determinante, 
die sich aus den Randbedingungen ergibt. Bei ein- 
gespannter Platte ist am Rande w — 0 und Be =!(. 
Auf diese Werte hat J keinen Einfluß, kann also bei 
der Bestimmung der Beullast bei beliebiger Beran- 
dung fortgelassen werden, nicht natürlich, wenn man 
die Auflagerreaktionen der ausgebeulten Platte be- 
rechnen will. Bei drehbar gelagerter Platte muß 
w=0( und M(s)—=0 sein, also ausführlicher ge- 
schrieben 


ds? 


wo das Glied mit «. als Faktor von J herstammt. 
Ist nun die Begrenzung aus Geradenstücken zu- 


sammengesetzt, so ist ar 0, außerdem ist für jede 


Also 


kann auch bei Bestimmung der Beullast der gerad- 
linig begrenzten Platte J im Ausdruck für die Form- 
änderungsarbeit fortgelassen werden, nicht aber bei 
der drehbar gelagerten, krummlinig begrenzten 
Platte. Während also die Beullast bei einer ebenen. 
beliebie geformten, eingespannten Platte und bei 
einer geradlinig begrenzten, drehbar gelagerten 
Platte von der Querkontraktionszahl « nur insofern 
abhängt. als diesesin dem Plattenfaktor N enthalten 
ist, hängt diese Beullast bei der nicht geradlinig be- 
grenzten, drehbar gelagerten, ebenen Platte außer- 
dem noch explizite von « ab. Das ist übrigens 
stets der Fall, wenn es sich um Platten handelt, 
deren Dicke veränderlich ist, etwa um abgestufte 
Platten, da dann J auf die Übergangsbedingungen 
an den Sprungstellen von Einfluß ist. 

Außerdem geht aus diesen Überlegungen hervor, 
daß man zur Bestimmung der Beullast einer nie ht 
geradlinig begrenzten, drehbar gelagerten Platte, 
diese nicht etwa ersetzen darf durch eine Platte, 
die durch ein passend gewähltes Polygon begrenzt 
ist, das sich aus Geradenstücken zusammensetzt. 
Das Gleiche gilt für die Berechnung der Durch- 
biegung krummlinig begrenzter, drehbar gelagerter 
ebener Platten, wie man aus der Gleichung Ia des 
Punktes 2 erkennt. Diese Tatsachen muß man bei 
einem Versuch, das Differenzenverfahren auf solche 
Platten anzuwenden, beachten. 


Dresden. 


2 
ebene gelagerte Platte am Rande = 


Willers. 962 


Bemerkungen zur graphischen iInte- 
gration. Bei der graphischen Integration einer 
gewöhnlichen Differentialgleichung 4’=1/(z, y) 
werden üblicherweise zwei Verfahren benutzt, die 
durch unmittelbare Übertragung der Formeln 2. Ord- 
nung von Runge, Heun und Kutta!) ins 
Zeichnerische gewonnen und meist in naheliegender 
Weise als Tangenten- bzw. Sehnenverfahren 
(im folgenden TV bzw. SV) bezeichnet werden. 
Beim TV gelangt man vom mg wei P (z,, Yo) 
mit der Schrittlänge k über den Zwischenpunkt 


zum Punkte Q*(z,+h, 
Y + %*) mit 


k* ist damit der zum TV gehörige Näherungswert 
für den wahren Ordinatenzuwachs 4y = k längs der 


1) Vgl. z. B. G. Schulz: rmunmiese zur prak- 
tischen Berlin 1937, S. 


Integralkurve durch P. Beim SV kommt man von 
P aus zum Punkte Q** h, k**) mit dem 
Näherungswert 


für k. Entwickelt man k* und k** nach Potenzen 
von Ah. so erhält man unter Benutzung des Operators 


- 
und 


Da die Entwicklung des genauen Ördinatenzu- 
wachses k 


2 
(5) 


1 
tz R? yo + + 
lautet, erkennt man. daß k* und k** bis auf einen 
Fehler 0(A®) mit k übereinstimmen. 

In der Vermutung, daß dadurch die Genauigkeit 
erhöht werden kann, ist nun in einer früheren 
Arbeit?) beiläufig vorgeschlagen worden, das TV 
in der folgenden Weise durcheine anschließende 
Iteration zu ergänzen. Man faßt das TV als 
ihren ersten Schritt auf und geht dann in der Weise 
weiter, daß man auf der y-Parallelen (z,—+ A) durch 
ein von Z ausgehendes Strahlenbüschel eine Folge 
von Punkten Q, = Q*, Q,, SER FRE . derart er- 
zeugt, daß jeweils ZQ,„ die dem Punkte Qu-.ı im 
Richtungsfeld der Differentialgleichung zugeordnete 
Richtung besitzt. Sind die Ko- 
ordinaten von Q,„. so ist also 


Nach den bekannten Sätzen über die Lösung einer 
Gleichung (z) dureh die Iteration (2, -,) 
konvergiert dieses Verfahren gegen die Lösung k 
der Gleichung 


wenn nur in der ganzen in Frage kommenden Um- 
gebung von Q(z,+h, y+ A) 


ist. Das bedeutet aber im wesentlichen die ge- 
läufige und im allgemeinen als erfüllt anzusehende 
Forderung nach der Beschränktheit von B R 

Wie festgestellt werden konnte, 
findet sich der orschlag desselben Iterationsver- 
fahrens bereits in einer kurzen Bemerkung bei 
Willers®), und zwar mit den Worten: „Noch ge- 


?2) H. Heinrich: Die graphische Integration von 
Differentialgleichungen ... mit nomographischen Hilfs- 
mitteln. Mitt. techn. Inst. Staatl. Tung-Chi-Univ., Woosung 
(China), Bd. I, Heft 1, S. 9. 


3) Fr. A. Willers: Methoden der praktischen Analysis. 
Berlin 1928, S. 306. 
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nauer wird die ...* Da hier 
jedoch ebenso wie in der oben angeführten Arbeit?) 
auf Genauigkeitsbetrachtungen verzichtet worden 
ist, das Verfahren aber offenbar, wie sich gezeigt 
hat, sehr naheliegt, dürfte es von Wert sein, die 
entsprechenden Untersuchungen im folgenden kurz 
nachzuholen und durch einige Bemerkungen zu er- 
gänzen. 

Zunächst ergibt sich für k, das ja im gegebenen 
Falle von A abhängt, aus (7) die Potenzreihen- 
entwicklung 


1 1 

und es zeigt sich, daß auch k nur bis auf einen 
Fehler % (A?) mit k übereinstimmt. Wollte man da- 
her von einer verbesserten Genauigkeit gegenüber 
dem TV oder dem SV sprechen, so könnte dies nur 
mit dem Hinweis darauf geschehen, daß sich die 
Entwicklung von %k infolge der Berücksichtigung 
des Gliedes f, D / im A?-Glied offensichtlich in engerer 
Weise an die Entwicklung von k anpaßt als die 
von k* und k**. Die Anwendung des geschilderten 
iterativen Verfahrens wird sich daher nach allem 
nur dann lohnen, wenn sich für jeden beliebigen 
Punkt die Feldrichtung einfach und genau genug 
bestimmen läßt. Ich kann, was diesen Punkt be- 
trifft, auf meine Ausführungen a.a.0.?) und in einer 
späteren Arbeit?) verweisen. 

Folgendes mag noch bemerkt werden. Die Be- 
rücksichtigung von /, D/ im h®-Glied von k ermög- 
licht es, sofort eine einfache lineare Verbindung 
von k, k*, k** anzugeben, deren Entwicklung nach 
Potenzen" von h bis auf einen Fehler 0 (A*) mit k 
übereinstimmt, nämlich 

3 2 

Oder man verbindet den zu dem iterativen Nähe- 

rungswert Ak gehörigen Punkt Q mit P und_geht 

von P aus in der zum Mittelpunkt M von PQ zu- 

geordneten Richtung bis zum Schnittpunkt Q*+** 

mit der y-Parallelen (2,+ Ah) vor. Dann ist das 
zugehörige 


R) 


und 


unterscheidet sich in der kA-Entwicklung ebenfalls 
erst im h*-Glied von Ah. In beiden Fällen muß je- 
doch dahingestellt bleiben, ob die dabei notwendigen 
Konstruktionen einfach genug sind, um in der 
Praxis der graphischen Integration noch verwend- 
bar zu sein. 

Man kann an das Iterationsverfahren noch von 
einem etwas anderen Gesichtspunkt herangehen. 
hl stellen dazu folgende Überlegungen an. 

Q(z+h, %-+ k) sei wie bisher der Punkt bei 

+ h, der genau auf der durch P gehenden Inte- 
% liegt. Die in P und Q an gelegten 
Tangenten mögen sich in Z(z,+ #9 h,n) schneiden. 


4 H. Heinrich: Ein an Ver- 
fahren zur Demusee von Richtungsfeldern ath. 
Bd. 3 (1938), S. 353 bis 389. 


(Dabei können wir ohne wesentliche Beschränkung 
der Allgemeinheit annehmen, daß 0O<9<1,d.h. 
daß Z „Zwischen“-Punkt ist. Das ist sicher immer 
dann der Fall, wenn die Krümmung von {% auf dem 
Weg von P nach @ ihr Vorzeichen nicht wechselt. 
Durch eine geeignete Schritteinteilung wird sich 
das stets erreichen lassen.) Für k ergibt sich dann 
die Gleichung 


k=dhh (13). 
Setzt man- zu ihrer Lösung die Iteration 
AÄn-ı). (14) 


an, so entspricht das sich daraus ergebende gra- 
phische Verfahren genau dem oben angegebenen 
mit dem einzigen Unterschied, daß der Ausgangs- 
unkt der Richtstrahlen nicht Z, sondern Z ist. 
ie Konvergenz ist gesichert, sobald in der ganzen 
in Frage kommenden Umgebung von Q 


ist, was im wesentlichen wieder auf die Beschränkt- 


heit von & hinauskommt. Wenn es daher ein — 


vielleicht auch iteratives — Verfahren gäbe, um 
den im allgemeinen nicht bekannten Zwischenpunkt Z 
zu bestimmen, so könnte man von ihm aus durch 
das angegebene Iterationsverfahren genau zu dem 
Bee Punkt Q der Integralkurve gelangen. 

olange jedoch Z nicht bekannt oder bestimmbar 
ist, ist man wieder auf Näherungen angewiesen. 
Aus (13) ergibt sich 


d ist mit k eine Funktion von A, für die wir die 


Reihe 


(16). 


(17) 


ansetzen können. Ist in etwas allgemeinerer Form 
als in (5) 


k=hy+ 
(18), 


wobei ,() die Ableitung niedrigster Ordnung > 2 
ist, die an der Stelle (z,, 4.) nicht verschwindet, 
so erhält man aus (16) nach einiger Zwischen- 
rechnung die für alle n>0 gültige Formel 


n+l1 
+ ("7") Ban . (19). 


Sie liefert für n=0 und n=1 


1) 


und daraus rekursiv d,, -.. 
wicklung von # beginnt daher mit 


h y“+D 


Die Reihenent- 


(20). 


Für den weitaus wichtigsten Fall, daß bereits 


Yo ist («=2), ergibt sich 


12 Y ” ® (20,1) 
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und, wenn P ein „gewöhnlicher“ Wendepunkt der 
Integralkurve mit 4," % 0(x=3) ist, 


3 
Es läßt sich leicht zeigen, daß das Arbeiten mit 
einem Zwischenpunkt, dessen ® mit einem Fehler 
0(h”) behaftet ist, im Ordinatenzuwachs k einen 
Fehler O(A”+*) nach sich zieht. Im Fall «=2 
müßte man daher, um den Fehler in k wenigstens 
auf O(A*) zu drücken, einen #-Wert benutzen, 
dessen A-Entwicklung sich frühestens im A°-Glied 
von (20,2) unterscheidet. Konstruktionen für solche 
Zwischenpunkte lassen sich zwar angeben, doch 
kann hier auf ihre Wiedergabe verzichtet werden, 
da sie nicht einfach genug erscheinen, um für die 
praktische Verwendung empfohlen werden zukönnen. 
Anschaulich gewinnt man ein Maß für dasin (20,1) 


. (20,2). 


wenn man den engen Zu- 


sammenhang zwischen %” und der Krümmung der 
Kurve und die Beziehung Er u ir | beachtet. 
Man kann auch davon ausgehen, daß der genaue 
Wert von # für die Integralkurve übereinstimmt 
mit dem entsprechenden für die eindeutig bestimmte 
Parabel 3. Ordnung, die durch P und Q hindurch- 
geht und dort mit der Integralkurve gemeinsame 
Tangenten besitzt. Ersetzt man diese Parabel 
näherungsweise durch diejenige, die durch P geht 
und dort mit der Integralkurve gemeinsame Ab- 
leitungen bis zur 3. Ordnung einschl. hat, und führt 
die entsprechende Rechnung durch, so erhält man 
als Näherungswert für # den Ausdruck 


% + hys 


nach Potenzen von h ebenfalls 
mit 5 13 + beginnt. 

Wichtig — weil in den bisherigen Darstellungen 
m. E. nicht genügend hervorgehoben — scheint mir 
ein besonderer Hinweis darauf zu sein, daß im 
Fall «=3, in dem P Wendepunkt der Integral- 
kurve ist, mit dem Näherungswert 3 für 9 ge- 
arbeitet werden muß und daß man dann — auch 
schon beim TV — sofort eine Übereinstimmung in 
k bis zum h®-Glied bekommt. 

Schließlich mag noch darauf aufmerksam gemacht 
werden, daß sich das Iterationsverfahren in der 
gleichen Weise wie bisher auch dann noch durch- 
führen läßt, wenn die Unterteilung des Richtungs- 
feldes nicht durch -Parallelen, sondern durch 
irgendeine andere Schar sich nicht überschneidender 
Kurven geliefert wird, wie es z. B. bei Verwendung 


von Leitkurven die Isopunktalen®) sind. Zur Ab- 
leitung der Konvergenzbedingeung kann man sich 
dabei mit Hilfe der vektorischen Darstellungsweise 
vom Koordinatensystem unabhängig machen. Die 
Rechnung selbst soll hier übergangen werden. 
Man gelangt zu folgendem Ergebnis. Sind p,q,.3 
die Ortsvektoren der Punkte P,Q,Z und w der 
Einheitsvektor der Feldrichtung. so daß mit dem 
Ortsvektor x des variablen Punktes X und der 
Bogenlänge s der Integralkurve als Parameter die 
Differentialgleichung in der Form HET er- 
scheint; ist weiter 8 die zu der Schar gehörige 
Kurve, auf der @ liegt, und ® die Bogenlänge 
längs $. so konvergiert das von Z aus angesetzte 
Iterationsverfahren dann, wenn für jeden Punkt 
X(g) auf 8 in der ganzen in Frage kommenden 
Umgebung von Q 


\d m! 
ZX 


ist. Das bedeutet. daß die Konvergenz um so 
besser ist, 


zg<i (21) 


a) je enger die Kurven & liegen (ZX klein); 
b) je weniger sich die Feldrichtung längs 8 


ändert ( om 


e) je schärfer die Feldrichtung die Kurve & 
schneidet (Winkel zwischen dem Tangenten- 


vektor Ss von 8 und dem Feldrichtungs- 


vektor nahe an 


Ist man mangels Bestimmbarkeit der genauen 
Lage von Z gezwungen, das_Iterationsverfahren 
näherungsweise von Punkten Z aus durchzuführen, 
so lassen sich auf diesen Fall die oben für die 
Unterteilung des Richtungsfeldes durch y- Parallelen 
angestellten Genauigkeitsbetrachtungen und ihre 
Ergebnisse weitgehend übertragen. 

Sind die Kurven & die Isoklinen des Richtungs- 
feldes, so reduziert sich das Iterationsverfahren 
auf dasTV. Sind sie außerdem zur y-Achse parallel, 
so handelt es sich um die Quadratur y=|f(z)dz, 
bei deren graphischer Durchführung man in dem 
bekannten Verfahren eine Möglichkeit zur Be- 
stimmung des Zwischenpunktes besitzt. 


Helmut Heinrich. 


| 
klein ); 


Breslau. 107 


5) Vgl. z. B. 
kurven in Richtungsfeldern. 
Bd. 19 (1939), S. 55. 


H. Heinrich: Allgemeines über Leit- 
Z. angew. Math. Mech. 
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C. B. BIEZENO, Prof. a. d. Techn. Hochschule 
Delft, R. GRAMMEL, Prof. a. d. Techn. Hoch- 
schule Stuttgart, Technische Dynamik, 
XII + 1056 S. m. 667 Abb. u. 5 Anhängen. Berlin 
1939, Verlag Julius Springer. Preis geb. 78 M. 

Die Verfasser, die einen großen Teil ihrer bis- 
herigen wissenschaftlichen Arbeit den Fragen der 
Elastizitätstheorie und der Dynamik, insbesondere 
ihren Anwendungen auf den Maschinenbau, gewidmet 
haben, behandeln in dem vorliegenden umfang- 
reichen Werk zahlreiche Probleme, die für den 
wissenschaftlich arbeitenden Ingenieur von Be- 
deutung sind. Der Begriff Dynamik wird hier 


weiter gefaßt als üblich — die Verfasser verstehen 
darunter ganz allgemein die Lehre von den 
Kräften —, so daß man in dem Buch vieles findet, 
was man dem Titel nach nicht erwartet. Nicht nur 
kinetische, sondern auch zahlreiche statische Pro- 
bleme, wie sie die technische Praxis der Mechanik 
stellt, werden darin behandelt. Natürlich mußte 
selbst für ein so umfangreiches Werk eine Aus- 
wahl getroffen werden. Daher haben sich die Ver- 
fasser in der Hauptsache auf Fragen beschränkt, 
die mit dem Maschinenbau, insbesondere mit dem 
Kraftmaschinenbau zusammenhängen. Vor allem 
sind ausführlich die Probleme behandelt, zu deren 
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Lösung einer der beiden Verfasser oder gar beide 
wesentliche Beiträge geliefert haben. Knapper be- 
trachtet sind Gegenstände, von denen bereits aus- 
führliche Darstellungen vorliegen, wie die Schalen- 
theorie, die Regelung usw. Es gibt aber nicht viel 
Probleme, die für den Ingenieur in den betreffenden 
Gebieten von Interesse sind, die man hier nicht 
findet. 

Vorausgesetzt wird an Mathematik und Mechanik 
im allgemeinen das, was etwa in den Kursus- 
vorlesungen der Technischen Hochschulen gebracht 
wird, nur hier und da einiges darüber hinaus, wie 
gewisse Kenntnisse aus der Theorie der Integral- 
gleichungen, der Funktionstheorie usw. Die Dar- 
stellung bevorzugt anschauliche Ableitungen 
und kommt auch sonst der Denkweise des In- 
genieurs entgegen, z. B. dadurch, daß — vielleicht 
abgesehen von der Inversion — keine konformen 
Abbildungen benutzt werden, sondern anstatt 
dessen mit angepaßten Koordinaten gearbeitet 
wird, ferner, daß nicht auf rein mathematisch 
interessierende Fragen, wie Konvergenzunter- 
suchungen, Existenzbeweise usw., eingegangen wird. 
Ausführlich wird natürlich die für alle Ausweich- 

robleme wichtige Eindeutigkeitsfrage der Lösungen 
behandelt. Im Vordergrund der Darstellung stehen 
stets die Methoden der Lösung, aber die Lösungen 
selbst sind auch immer so weit geführt, daß zahlen- 
mäßige Resultate aus den gefundenen Formeln mit 
erträglichem Rechenaufwand abgeleitet werden 
können. Oft sind rechnerisch bequeme Näherungs- 
formeln, oft auch Rechenschemata für die Durch- 
führung der Zahlenrechnungen gegeben. Da für 
manche Leser nur einzelne Abschnitte von Interesse 
sein werden, haben die Verfasser die Darstellung 
so gestaltet, daß ein Leser, der schon etwas ver- 
traut mit dem Gegenstand ist, sich in die einzelnen 
Kapitel einarbeiten kann, ohne von den vorher- 
gehenden genauer Kenntnis genommen zu haben. 
Vielfach wird das Buch auch als Nachschlagewerk 
benutzt werden. Das ist auf Grund der Dar- 
stellungsart und des ausführlichen Inhaltsverzeich- 
nisses möglich, doch würde die Brauchbarkeit für 
diesen Zweck durch ein leider fehlendes, eingehen- 
des Sachverzeichnis erhöht werden. 

Den reichen Inhalt im einzelnen aufzuführen, 
erlaubt leider der hier zur Verfügung stehende 
Platz nicht; nur ein kurzer Überblick sei gegeben. 
Der Stoff ist in vier Abschnitte gegliedert. Im 
ersten werden die Grundgesetze und die allgemeine 
Theorie der Elastomechanik behandelt. Hier findet 
man unter anderem den klassischen Gleichungen 
des elastischen Gleichgewichtes einen Satz von 
Gleichungen an die Seite gestellt, die den neu- 
tralen Gleichgewichtszustand, also die Grenzlage, 
bei der der Gleichgewichtszustand bei weiterem 
Wachsen der Belastung aufhört eindeutig zu sein, 
allgemein erfassen. Die hier gewonnenen Gleichun- 
gen finden zwar bei den späteren Untersuchungen 
technischer Ausweichprobleme keine Verwendung, 
doch werden sie als Grundlage für eine verfeinerte 
Behandlung gute Dienste leisten können. Ferner 
werden hier ausführlich die verschiedenen Lösungs- 
methoden, insbesondere auch die Näherungsverfahren 
für Rand- und Eigenwertprobleme, die experimen- 
tellen Methoden zur Spannungsermittlung, wie die 
verschiedenen Analogiemethoden dargestellt. Im 
zweiten Abschnitt werden einzelne Maschinenteile 
wie Stab und Welle, Feder und Ring, Platte und 
Schale, ferner in einem besonderen Kapitel die Aus- 
weichprobleme, also die Knickprobleme, das Durch- 
schlagen gekrümmter Stäbe und Platten und das 
Umkippen von Ringen untersucht. Der dritte Ab- 
schnitt wendet sich Fragen zu, die mit dem Dampf- 
turbinenbau zusammenhängen: rotierende Scheiben, 
Dehnung, Biegung und Schwingungen der Schau- 
feln, kritische Drehzahlen der einfach und mehrfach 
besetzten Wellen usw. Der letzte Abschnitt endlich 


beschäftigt sich mit Problemen, die bei Bau und 
Untersuchung von Brennkraftmaschinen auftreten, 
so mit dem Massenausgleich bei den verschiedenen 
Motoren, dem Leistungsausgleich bei Ein- und Mehr- 
zylindermaschinen, den rehschwingungen ein- 
facher und verschieden zusammengeschalteter Ma- 
schinen, ferner mit den erzwungenen Schwingungen, 
der Resonanz und Scheinresonanz, dem Einfluß der 
Nebentorsion und dem veränderlicher Drehmassen. 

Das wertvolle Buch, das den umfangreichen Stoff 
in sehr geschickter Anordnung und Darstellung be- 
handelt, wird auf viele Jahre hinaus nicht nur allen, 
die auf diesen Gebieten arbeiten, ein zuverlässiger 
Ratgeber und ein unentbehrliches Nachschlagewerk 
sein, sondern auch vielen als ausgezeichnete Ein- 
führung in die verschiedenen Gebiete der tech- 
nischen Dynamik dienen. 


Dresden. Willers. 133 


Dr. L. FÖPPL, o. Prof. a. d. Techn. Hochschule 
München, Aufgaben austechnischer Me- 
chanik, Unterstufe, Statik, Festig- 
keitslehre, Dynamik. VII+202 S. m. 
317 Abb. München u. Berlin 1939, Verlag R. Olden- 
bourg. Preis kart. 10 M. 


Das Buch enthält zahlreiche, instruktive Aufgaben 
zur Statik, Festigkeitslehre und Dynamik, die aus 
der Ingenieurpraxis geschöpft und für den Unter- 
richt an der Technischen Hochschule aufbereitet 
wurden, sowie ihre Lösungen. Sicherlich wird auch 
die vorliegende zweite Auflage, die sich von der 
ersten, 1929 erschienenen nur wenig unterscheidet, 
gern von Lernenden und Lehrenden zur Hand ge- 
nommen werden. Man vergleiche die Besprechung 
Be Seen Auflage in Band X der ZAMM (1980), 

Freiberg (Sa.). Grüß. 134 

Dr.-Ing. RUDOLF JAESCHKE, Flugzeug- 
bereehnung. Bd.I, Grundlagen der Strömungs- 
lehre und Flugmechanik. 2. Aufl. 174 S. m. 88 Abb. 
u. 21 Zahlentaf. Bd. II, Bearbeitung von Entwürfen 
und Unterlagen für den Festigkeitsnachweis. 2028. 
m. 64 Abb. u. 38 Zahlentaf. München u. Berlin 1939, 
u R. Oldenbourg. Preis geb. in einem Bd. 


Wie die bekannte Misessche Fluglehre wendet 
sich das vorliegende Buch auch „an alle diejenigen, 
die sich aus beruflichen oder persönlichen Inter- 
essen mit Fragen des Flugzeugbaus beschäftigen 
wollen“. Trotzdem ist das Buch in erster Linie ge- 
schrieben für Studierend e des Flugzeugbaus: 
Es ist hervorgegangen aus Kursen des Ver- 
fassers an einer Ingenieurschule. Gegliedert ist es 
in zwei Bände, die zusammen 6 Hauptteile um- 
fassen: Strömungslehre, Flugleistungsberechnung, 
Flugzeugstabilität, Bauarten, Entwerfen von Flug- 
zeugen, Belastungsannahmen. Wie das nicht 
anders sein kann, deckt sich der Stoff zum großen 
Teil mit dem in der „Fluglehre“ behandelten — 
nur liegt, dem Titel des Buches entsprechend, das 
Hauptgewicht hier sehr viel mehr auf den Berech- 
nungsbeispielen. So sind in der „Fluglehre“ alle 
Flugzeugteile in etwa gleicher Ausführlichkeit und 
immer unter Hervorhebung des physikalisch Wich- 
tigen behandelt, während im vorliegenden Werk 
z. B. auf die Flugmotoren nur soweit eingegangen 
wird, als dies für die Beschaffung von Entwurfs- 
unterlagen wünschenswert erscheint. Nur die 
Strömungslehre ist relativ breit abgehandelt, aller- 
dings ist gerade in diesem Abschnitt die Ausdrucks- 
weise gelegentlich unscharf: So wenn nach S. 21 
die Zentrifugalkräfte der abgelenkten Luftteilchen 
den aus der Bernoullischen Gleichung folgenden 
Auftrieb „unterstützen“ (die B.sche Gleichung gibt 
den Auftrieb ja in aller Vollständigkeit wieder) 
oder wenn es auf S. 31 heißt „Die Zirkulation reißt 
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vollständig ab“, oder wenn in Abb. 20 die Prandtl- 
sche Kontrollinie nach Art der Stromlinien einen 
Pfeil trägt. Störender noch als diese gelegent- 
lichen Flüchtigkeiten scheint dem Referenten die 
Aufnahme des (in sich natürlich berechtigten) 
letzten Abschnittes: „Berechnung der Luftkräfte“ in 
das Buch; denn dessen mathematische Anforderun- 
gen gehen weit hinaus über die im übrigen Buch 
gestellten (s. z. B. die Einführung des Begriffs der 
Beschleunigung S. 10), so daß die Leser, für die 
das übrige Buch doch bestimmt ist, den Entwick- 
lungen dieses Abschnittes wohl kaum werden folgen 
können. — Trotz dieser Einwände im Ganzen ein 
Buch, das durch seine übersichtliche Gliederung und 
das darin gesammelte Zahlenmaterial dem Lernen- 
den auf dem Gebiete des Flugzeugentwurfs gute 
Dienste zu leisten vermag. 


Berlin-Adlershof. Marguerre. 146 


Dr.-Ing. F. RÖTSCHER, Prof. a. d. Techn. Hoch- 
schule Aachen, und Dr.-Ing. R. JASCHKE, Assist. 
a. d. Techn. Hochschule Aachen, Dehnungs- 
messungen und ihre Auswertung. 
VI+ 121 S. m. 191 Abb. und 1 Tafel. Berlin 1939, 
Verlag Julius Springer. Preis brosch. 16,80 M. 


Die durch Rechnung gewonnenen Ergebnisse der 
Ermittlung der Spannungsverteilung sind oft — aus 
verschiedenen Gründen — nicht ausreichend, um 
die einwandfreie Bemessung technischer Konstruk- 
tionen hinsichtlich der Bedingungen des gesicher- 
ten Materialverhaltens zu gewährleisten. Die 
Rechnung muß dann durch Messungen ergänzt 
werden. Da aber die Spannungen nicht unmittel- 
bar, sondern nur auf dem Wege über die Dehnun- 
gen ermittelt werden können, wird hier einerseits 
eine übersichtliche Beschreibung einer Anzahl von 
ausgeführten und bewährten Dehnungsmessern ge- 

eben, andrerseits aber — mit großer Ausführlich- 
Keit — die Frage behandelt, wie aus Dehnungs- 
messungen einwandfreie Rückschlüsse auf die Span- 
nungsverteilung gewonnen werden können. Dabei 
werden in umfassender Weise die Mohrschen Kreise 
herangezogen, die sich auch auf anderen Gebieten 
der technischen Mechanik als überlegenes Darstel- 
lungsmittel für Tensoren zweiter Stufe in zwei und 
drei Dimensionen bewährt haben. — Die Behand- 
lung von — anderswo selten gebrachten — Sonder- 
fragen, wie die Zusammensetzung von Verformungs- 
und Spannungszuständen, die der Fehler und Stö- 
rungen bei Dehnungsmessungen u. v. a., sowie auch 
die vielen lehrreichen Abbildungen erhöhen noch 
den Wert der sehr anschaulich geschriebenen und 
durch seine besondere Fragestellung und Darstel- 
lungsweise außerordentlich aufschlußreichen Büch- 
leins. 


Karlsruhe. Th. Pöschil. 152 


Dr. GUSTAV MESMER, Spannungsoptik. 
XI + 222 S. m. 197 zum Teil farbigen Abb. Ber- 
lin 1939, Verlag Julius Springer. Preis geb. 30 M. 

Mit Rücksicht auf die zunehmende Bedeutung der 
Spannungsoptik für die Praxis ist das vorliegende 
Buch sehr zu begrüßen. Es stellt gegenüber dem 
in deutscher Sprache bisher allein vorhandenen 
Lehrbuch über Spannungsoptik v. L. Föppl und H. 
Neuber eine wertvolle Ergänzung und Bereicherung 
dar. Besonders hervorzuheben sind hierbei die zahl- 
reichen gut ausgeführten praktischen Beispiele; 
ferner die Zusammenstellung der verschiedenen Ge- 
räte und der meisten Lösungsverfahren von prak- 
tischer Bedeutung. Auch der umfangreiche Litera- 
turnachweis am Schluß des Buches stellt eine wert- 
volle Zusammenstellung für den Forscher dar. Das 
erste Kapitel ist der üblichen Darstellung des 
ebenen Spannungszustandes, das zweite den 
optischen Grundlagen der Spannungsoptik gewid- 


met. Es folgen im 3. und 4. Kapitel die Lösungs- 
verfahren der Spannungsoptik, deren Bedeutung in 


re en Hinsicht gut abgeschätzt worden sind. 
as 5. Kapitel, das sich mit den erforderlichen Ge- 
räten beschäftigt, sowie das 6. Kapitel, das die Ver- 
suchstechnik im einzelnen behandelt, enthält 
namentlich für solche, die zu Forschungszwecken 
oder für die Praxis ein spannungsoptisches Labo- 
ratorium einrichten und mit solchen Versuchen be- 
ginnen wollen, wertvolle Hinweise. Das 7. Kapitel 
geht auf Sonderfragen der Spannungsoptik ein, die 
erst in den letzten Jahren behandelt worden sind, 
nämlich auf dynamische Spannungsaufgaben und 
ferner auf das sehr wichtige und zukunftsreiche Ge- 
biet der räumlichen Spannungsoptik. Das Buch 
schließt mit dem 8. Kapitel, das eine reiche Über- 
sicht über verschiedene auf spannungsoptischem 
Wege gelöste praktische Aufgaben gibt. — 

Das Buch kann jedem, der sich mit Spannungs- 
optik beschäftigen will, wärmstens empfohlen 
werden. 


München. L. Föppl. 144 


Dr. A. THUM, Prof. a. d. Techn. Hochschule 
Darmstadt, und Dr.-Ing. K. FEDERN, Span- 
nungszustand und Bruchausbildung., 
anschauliche Darstellung der span- 
nungsmechanischen Grundlagen der 
Gestaltsfestigkeit und der Gesetz- 
mäßigkeiten der Bruchausbildung. 
V+78S. m. 8 Abb. Berlin 1939, Verlag Julius 
Springer. Preis brosch. 9,60 M. 


Auf Grund der Ergebnisse zahlreicher Dauer- 
versuche, die von den Verfassern an der Staat- 
lichen Materialprüfungsanstalt Darmstadt durchge- 
führt wurden, werden im vorliegenden Bande für die 
Gesetzmäßigkeiten der Bruchausbildung und der 
Bruchgefahr spannungsmechanische Deutungen ge- 
geben. Dabei knüpfen die Verfasser an neuere 
Erkenntnisse der theoretischen Festigkeitsforschung 
(vgl. H. Neuber, Kerbspannungslehre, Berlin 1937, 
Julius Springer) an und gelangen in einer Reihe 
von Fällen zu wichtigen Folgerungen. Das unter- 
schiedliche Verhalten der Werkstoffe gegenüber 
„Gewaltbruch“, „Zeitbruch“ und „Dauerbruch“ wird 
an Hand von vielen Bruch- und Gefügebildern 
unter Berücksichtigung von Beanspruchungsart und 
Spannungszustand eingehend erläutert. Der Ein- 
fluß des Spannungszustandes auf den Bruchvor- 
gang, der an Hand der vielen Versuchsergebnisse 
besonders eindringlich hervortritt, ist ein Beweis 
der großen Bedeutung rechnerisch erfaßbarer Span- 
nungszustände für die Dauerbruchforschung. 


Braunschweig. H. Neuber. 145 


Oberbaurat Dr.-Ing. ehr. FRITZ von EMPERGER, 
Stahlbeton mit vorgespannten Zu- 
lagen aus höherwertigem Stahl. VI+ 
38 S. m. 22 Abb. Berlin 1939, Verlax Wilhelm 
Ernst & Sohn. Preis brosch. 4,40 M. 


Während die elementare Festigkeitslehre für 
Verbund-Konstruktionen aus Beton und Stahl von 
vornherein mit der gerissenen Zugzone des Quer- 
schnitts rechnet, beschäftigt sich der Verfasser in 
seiner Arbeit mit Maßnahmen zur Vermeidung von 
Rissen und zur Erhöhung der Wirtschaftlichkeit 
dieser Baukonstruktionen durch die Verwendung 
von Stahl von höherer Festigkeit und durch die 
Vorspannung der Rundeisenbewehrung. Er er- 
örtert in Verbindung damit Maßnahmen, dieses Ziel 
namentlich durch die gleichzeitige Verwendung ver- 
schiedener Stahlsorten zu erreichen, und begründet 
seine Vorschläge durch die Ergebnisse von Ver- 
suchen, die insbesondere der Haftfestigkeit von 
dünnen, vorgespannten Drähten im Beton gelten. 
Diese Untersuchungen bilden einen wertvollen 
empirischen Beitrag zur Förderung der gegenwärti- 
gen Bestrebungen im Eisenbetonbau. 


Dresden. K. Beyer. 141 
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Dr.-Ing. FRITZ LEONHARDT, Anleitung für 
die vereinfachte Trägerrostberech- 
nung mit Hilfstafeln, Formeln und 
Beispielen. III + 104 S. m. 86 Abb. Berlin 
age Wilhelm Ernst & Sohn. Preis brosch. 
7,50 M. 


Der Verfasser behandelt die Entwicklung von 
Modellen von Trägerrosten und die Ausmessung 
des Verschiebungszustandes unter ausgezeichneten 
Belastungseinheiten, um daraus die Einflußlinien 
der Verbindungskräfte an den Knotenpunkten des 
Rostes nach dem Maxwellschen Gesetz abzuleiten. 
Das Verfahren, das bereits anderweit ausgeprobt 
ist, eignet sich hier insbesondere zur statischen 
Untersuchung von Trägerrosten mit unregelmäßiger 
oder schiefer Zuordnung der Träger, da sich die 
analytischen Untersuchungen dann nicht durch 
Symmetriebetrachtungen und Belastungsumord- 
nung vereinfachen lassen. In dem 2. Teil des 
Buches werden die Verbindungskräfte an den 
Knoten regelmäßiger Trägerlagen mit einem Quer- 
träger angegeben und für mehrere Querträger ge- 
schätzt. Diese Ergebnisse sind zwar bereits be- 
kannt, aber hier übersichtlich für den Handgebrauch 
am Konstruktionstisch bearbeitet, zusammengestellt 
und durch Beispiele erläutert worden. Ein Schluß- 
abschnitt enthält wertvolle Richtlinien für die Aus- 
bildung von Trägerrosten in Stahl- und Eisenbeton. 


Dresden. K.Beyer. 142 


Dr.-Ing. HEINRICH PRESS, Der Boden als 
Baugrund. (Mitteilungen aus dem Gebiete des 
Wasserbaues und der Baugrundforschung, Heft 9.) 
V-+49S.m. 64 Abb. Berlin 1939, Verlag Wilhelm 
Ernst & Sohn. Preis br. 5,80 M. 


Das Mitteilungsheft aus dem Gebiete der Bau- 
grundforschung bringt nicht, wie man zunächst er- 
wartet, neue Forschungsergebnisse, sondern faßt 
die in den letzten 10 Jahren von Preß in ver- 
schiedenen Zeitschriften veröffentlichten Versuchs- 
ergebnisse und sonstigen Abhandlungen zusammen. 
Da diese nur einige kleine Teilgebiete des großen 
Themas „Der Boden als Baugrund“ betreffen, hat 
der Verfasser die Lücken durch Auszüge aus 
Normen, Richtlinien und Baugrundlehrbüchern er- 
ea Da er nun aber nur sich selbst zitiert, 
könnte ein unbefangener Leser den Eindruck ge- 
winnen, als sei die neuere Bodenmechanik von 
Preß aufgebaut worden. Demgegenüber muß fest- 
gehalten werden, daß selbst die eigenen Preßschen 
Versuche meist nur eine Wiederholung älterer Ver- 
suche darstellen, also meist ohne Originalität oder 
Priorität sind. 

Auf den ersten 10 Seiten behandelt das Heft die 
Bodenbenennung (Din 4022), allgemeine, 
geophysikalische und bautechnische Bodenauf- 
schlüsse und die Bodenentnahme (Din 


4021). Die restlichen 40 Seiten sind der Boden- 
tragfähigkeit gewidmet. Hiervon stellen 


wiederum 10 Seiten zunächst einen kurzen lehr- 
buchmäßigen Auszug aus der Bodenphysik dar, 
dann folgen weitere 10 Seiten Versuchsergebnisse 
von Probebelastungen des Verfassers mit kleinen 
Lastflächen, denen sich Ergebnisse von Pfahlprobe- 
belastungen anschließen (S. 31—37). Es folgt ein 
kurzer lehrbuchmäßiger Auszug aus der Grund- 
wasserlehre. Die letzten 10 Seiten bringen die 
Druckverteilung im Baugrunde und die Setzungs- 
ermittlung. Bei der Druckverteilung werden in 
rein autodidaktischer Form auch nur eigene Ver- 
suche angeführt, der Rest ist wiederum ein Aus- 
zug aus der Setzungslehre ohne Anführung der be- 
nutzten Quellen. 

Obwohl Preß von dem verstorbenen Prof. Dr. 
Kögler schon einmal darauf aufmerksam gemacht 
worden ist, daß es doch bei wissenschaftlichen Ver- 
öffentlichungen üblich sei, auch auf fremde Arbeiten 


hinzuweisen, besonders wenn sie genau das 
gleiche behandeln, aber zeitlich früher ausgeführt 
worden sind, hat er wiederum an seiner ungewöhn- 
lichen einseitigen Darstellungsweise festgehalten. 

Preß fordert schließlich die Fachkreise eifrig zu 
Setzungsmessungen auf, veröffentlicht in seinem 
Büchlein aber nicht ein einziges eigenes Ergebnis, 
das doch überzeugender werben müßte als bloße 
Worte. 

Ohne den Wert der Preßschen Versuche, die er 
mit Hingabe meist zur Bestätigung früherer ähn- 
licher Versuche durchgeführt hat, irgendwie herab- 
setzen zu wollen, mußte doch die Art und Form 
der Darstellung im Interesse des unbefangenen 
Lesers festgestellt werden. Ob der in der Bau- 
praxis stehende Ingenieur den Überblick über das 
Thema gewinnt, den Preß als Ziel hinstellt, muß 
bei dem einseitigen Darstellungsverfahren füglich 
bezweifelt werden. 

Naumburg (Saale). Dr. Scheidig. 139 

Dr.-Ing. habil. KURT RAUH, a. o. Prof. a. d. 
Techn. Hochschule Aachen, Praktische Ge- 
triebelehre. II. Bd. VII +298 S. m. 709 
Abb. Berlin 1939, Verlag Julius Springer. Preis 
geb. 29,40 M. 

Wenn in dem Gelenkviereck einzelne Gelenke in 
Geradführungen umgebildet werden, so können dar- 
aus zahlreiche andere Getriebe abgeleitet werden. 
Geschieht dies mit drei Gelenken, so entsteht die 
dreigliedrige Keilkette, von der durch Spezialisie- 
rung drei weitere Getriebegruppen abgeleitet wer- 
den können: die Schubgetriebe, Sperrungen oder 
Kupplungen und die Führungsgetriebe, die in die- 
sem Bande eingehend behandelt werden. Die bis- 
her benutzte Reuleauxsche Einteilung der Getriebe 
erweist sich dabei als unzulänglich und ergibt nur 
einige der möglichen Stufen eines systematischen 
Getriebeaufbaues. Zahlreiche Einzelfragen der prak- 
tischen Getriebelehre, die sich auf den Entwurf und 
die einwandfreie Gestaltung und Herstellung ge- 
räuscharmer Getriebe, auf die Kupplungen und 
stufenlosen Schaltungen u. dgl. beziehen, werden 
besonders behandelt. Das Buch enthält zahlreiche, 
übersichtlich angeordnete Abbildungen und prak- 
tische Hinweise und dürfte seiner elementaren Be- 
schaffenheit halber, die mathematische Entwicklun- 
en fast völig vermeidet, dem praktischen Getriebe- 
onstrukteur willkommen sein. 


Karlsruhe. 


Dipl.-Ing. W. BÜTTNER, Die Entropie, 
physikalische Grundlagen und.tech- 
nische Anwendungen. 140 S. m. 86 Abb. 
1939, VDI-Verlag G.m.b.H. Preis brosch. 
9M. 


Th. Pösehl. 153 


Die Entropie ist eine Zustandsgröße, mit der der 
Ingenieur sehr viel zu tun hat: Sie ist ja die Ab- 
szisse der meisten Mollierdiagramme. Trotz dieses 
häufigen Gebrauches stehen sehr viele Ingenieure 
der Entropie mit einer gewissen Scheu gegenüber 
— vielleicht weil sie bei den üblichen Ableitungen 
dieses Begriffes die Gegenständlichkeit vermissen. 

Hier will die Arbeit von Büttner weiterhelfen. 
In dem etwa 30 Seiten umfassenden ersten Teil der 
Schrift wird die Zustandsgröße Entropie abgeleitet. 
Dabei wird auch auf den Zusammenhang zwischen 
Entropie und Wahrscheinlichkeit eingegangen. 

Der zweite umfangreichere Teil des Büchleins 
une in klarer Weise die Grundlagen und die An- 
wendung der Entropieschaubilder auf Wärmekraft- 
und -arbeitsmaschinen. 

Im Anhang wird noch das Entropiediagramm der 
festen Körper behandelt und gezeigt, wie es bei 
Festigkeitsrechnungen zur Bestimmung des Sicher- 
heitsgrades benutzt werden kann. 


Dresden. Mehlig. 135 
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WILHELM WEYDANZ VDI, Die Vorgänge 
in Strahlapparaten. Mitteilungen aus dem 
Maschinenlaboratorium der Technischen Hochschule 
Karlsruhe. (Beiheft zur Zeitschrift für die gesamte 
Kälte-Industrie, herausgegeben von Prof. Dr.-Ing. 
R. Plank, iteiie 2, Heft 8.) 139 S. m. 53 Abb. 
Berlin 1939, VDI-Verlag G.m.b.H. Preis brosch. 
3M. 


Es ist ein Zeichen für die Wichtigkeit, die den 
Strahlapparaten beigemessen wird, daß ganz kurz 
nach der Arbeit von Flügel über die Berechnung 
von Strahlapparaten eine umfangreichere Arbeit er- 
scheint, die das Problem von einer völlig anderen 
Seite anpackt. Es wird in der Arbeit von Weydanz 
der Versuch unternommen, Berechnungsunterlagen 
für den Entwurf von Strahlern auf Grund von Mes- 
sungen zu gewinnen, die zum Teil von Verfasser 
selbst durchgeführt wurden. Bei der Berechnung 
wird der Strömungsvorgang im Strahler aufge- 
spalten in drei Teile, nämlich die Strömung in der 
I>üse, im Mischraum und im Diffusor. Für diese 
drei Strömungen werden die Verluste und die sich 
daraus ergebenden Einzelwirkungsgrade berechnet. 
Am wertvollsten ist dabei die Untersuchung über 
den Wirkungsgrad der Düse auf Grund des be- 
kannten Ansatzes über die Wandreibung. Die Er- 
gebnisse dieser Rechnung haben auch dem Dampf- 
turbinenbauer mancherlei zu sagen. 

Das Ergebnis der Arbeit ist eine Berechnungs- 
anweisung, die es gestattet, beim Neuentwurf von 
Strahlern die Querschnitte zu bestimmen und die 
richtigen Druckverhältnisse zu wählen. 

Der Verfasser bezeichnet selbst die zahlenmäßi- 
gen Ergebnisse seiner Arbeit als eine erste An- 
näherung und die ganze Untersuchung als einen 
ersten grundsätzlichen Schritt zu allgemeinen Be- 
rechnungsmethoden für Strahlapparate. Selbst unter 
Berücksichtigung dieser Einschränkung müssen aber 
einige grundsätzliche Bedenken erhoben werden, 
die zur Vorsicht mahnen. 

1. Bei den Versuchen wurden nur die Drucke und 
Temperaturen vor und hinter dem Strahler sowie 
die durchfließenden Stoffmengen gemessen. Mes- 
sungen zur Klarstellung der Strömungsvorgänge 
innerhalb des Strahlers wurden überhaupt nicht ge- 
macht. Infolgedessen beruhen sämtliche Aussagen 
des Verfassers über die Strömungen im Strahler 
nur auf Annahmen. Bei der heutigen Versuchs- 
technik wäre es möglich gewesen, wertvolle Auf- 
schlüsse durch verhältnismäßig einfache Messungen 
zu gewinnen. 

2. Die Versuche wurden nur an fertigen Strahlern 
durchgeführt, an deren Abmessungen nichts ge- 
ändert werden konnte. Vor allem wurde nur die 
gerade vorhandene Form des Mischraums unter- 
sucht, trotzdem anzunehmen ist, daß diese auf 
Wirkungsgrad und Kennlinie des Strahlers einen 
sehr großen Einfluß hat. 

3. Die strömungstechnischen Erklärungen des 
Verfassers sind an mehr als einer Stelle anfechtbar 
und lassen sich nicht mit dem heutigen Stand der 
Strömungslehre in Einklang bringen. 

Aus diesem Grunde muß es selbst unter den 
obengenannten Einschränkungen als fraglich er- 
scheinen, ob die Arbeit den Ingenieur wirklich in 
der Vorausberechnung von Strahlern wesentlich 
weiterführt. 


Augsburg. E. Sörensen. 136 


ARNOLD SOMMERFELD, Prof. a. d. Universität 
München, Atombau und Spektrallinien. 
I. Bd. 2. umgearb. u. erweit. Aufl. d. „Wellen- 
mechanischen Ergänzungsbandes“. XI-+ 820 S. m. 
62 Abb. Braunschweig 1939, Verlag Vieweg & Sohn. 
Preis geb. 35 M. 


Diese langerwartete zweite Auflage wird von 
Physikern und Mathematikern begrüßt werden vor 


allem als eine einzigartige Darstellung der mathe- 
matischen Hilfsmittel der Wellenmechanik. Diese 
Hilfsmittel sind vollständig entwickelt und an der 
konkreten Aufgabe anknüpfend weitgehend allge- 
mein durchgeführt. Das Grundsätzliche in der 
physikalischen Auffassung und die Begründung der 
Ansätze durch allgemeine Erfahrungstatsachen ist 
dabei etwas in den Hintergrund gestellt. Dafür 
sind jedoch eine Reihe von Fragen, die für die Be- 
obachtung wichtig sind (Photo-Effekt, kontinuier- 
liches Röntgenspektrum, Compton - Effekt) voll- 
ständig durchgerechnet. Außer den genannten 
Fragen sind gegenüber der ersten Auflage wesent- 
lich ausführlicher dargestellt einige Eigenschaften 
des Wasserstoffspektrums (dabei auch die hypergeo- 
metrischen Funktionen) und vor allem die Diracsche 
Theorie des Elektrons unter Benutzung der Gruppe 
der hyperkomplexen Einheiten (statt der beson- 
deren Diracschen Matrizen). 


Leipzig. F. Hund. 122 


Ferner sind bei der Schriftleitung folgende 
Bücher ge (ausführliche Besprechung 
bleibt vorbehalten): 


C. CRANZ und O. von EBERHARD, Die neu- 
zeitliche Entwicklung der Schuß- 
waffen. 2. erweiterte Aufl. 72 S. m. 23 Bildern. 
Berlin 1939, VDI-Verlag G.m.b.H. Preis brosch. 
2M. 


Dr. se. techn. ALFRED MANGER, Der durch- 
laufende Balken auf elastisch dreh- 
baren und elastisch senkbaren Stützen 
einschließlich des Balkens auf stetiger 
elastischer Unterlage. (Mitteilungen aus 
dem Institut für Baustatik a. d. Eidg. Technischen 
Hochschule in Zürich, herausgegeben von den Pro- 
fessoren Dr. M. Ritter und Dr. F. Stüssi, Mit- 
teilung 10.) 170 S. m. 77 Fig., 25 Tab., 8 graph. 
Beilagen. Zürich und Leipzig 1939, Verlag A.-G. 
Gebr. Leemann & Co. Preis brosch. 7,50 M. 


II Vo Convegno Interprovineiale Di 
Aerotecnica, Neapel, 19. bis 22. Mai 1938. 
466 S. Rom 1939, Istituto Poligrafico Dello Stato. 


Dr. HERMANN DIESSELHORST, Prof. em. d. 
Physik a. d. Techn. Hochsch. in Braunschweig, 
Magnetische Felder und Kräfte mit 
einer Übersicht über die Vektoren- 
rechnung, eine mit Zusätzen versehene Sonder- 
ausgabe des Beitrages „Elektrodynamik“ aus Hand- 
buch der Elektrizität und des Magnetismus, Bd IV 
1920). VI + 216 S. m. 56 Abb. Leipzig 1939, Ver- 
lag Johann Ambrosius Barth. Preis kart. 8 M. 


Prof. Dr. F. RAUSENBERGER, Theorie der 
Rohrrücklaufgeschütze. 2. Aufl. bearb. v. 
Prof. Dr. O. von Eberhard. IX + 175 S. m. 64 Abb. 
Berlin 1939, Verlag Julius Springer. Preis geb. 
16,20 M. 


Professor Dipl.-Ing. RUÜDOLF VOGDT, Formel- 
sammlung zur Technischen Mechanik 
(Statik und Dynamik des Maschinen- 
baus). II. verb. Aufl. VIT+6 S. m. 132 Abb. 
Leipzig 1939, Dr. Max Jänecke Verlagsbuchhand- 
lung. Preis brosch. 0,00 M. 


Dr. K. SCHWIDEFSKY, wissenschaftl. Mitarbei- 
ter d. Opt. Werke Zeiß, Jena, Einführungin 
die Luft- und Erdbildmessung. 2. erw. 
u. verb. Aufl. m. 73 Abb., drei schwarzen und zwei 
farbigen Tafeln im Text, einer schwarzen Tafel, 
einer farbigen Brille und zwei Stereobildern im An- 
hang. IV + 137 S. Leipzig und Berlin 1939, Verlag 
B. © Teubner. Preis geb. 8 M. 
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Dr. phil. FRIEDRICH DRENCKHAKHN, Prof. a. d. 
Hochschule f. Lehrerbildung, Rostock, unter Mit- 
arbeit von Dr.-Ing. habil. ULRICH GRAF, Prof. a. 
d. Techn. Hochschule Danzig, Geländemathe- 
matik. VII+212 S. Langensalza 1939, Verlag 
Julius Beltz. Preis geb. 7 M. 


Dr. phil. HEINRICH HECHT, Schalt- 
schemata und Differentialgleichun- 
gen elektrischer und mechanischer 
Schwingungsgebilde VI+135 m. 
51 Abb. Leipzig 1939, Verlag Johann Ambrosius 
Barth. Preis brosch. 9,60 M. 


Professor Dr. phil. Dr.-Ing. e. h. KARL WILLY 
WAGNER, ÖOperatorenrechnung nebst 
Anwendungen in Physik und Technik, 
nach Vorträgen, veranstaltet durch den Bezirk Berlin 
des Verbandes Deutscher Elektrotechniker — vor- 
mals Elektrotechnischer Verein e. V. — in Ge- 
meinschaft mit dem Außeninstitut der Technischen 
Hochschule Berlin. XIV +448 S. m. 126 Abb. 
Leipzig 1940, Verlag Johann Ambrosius Barth. 
Preis brosch. 27,60 M. 


Dr.-Ing. FRITZ WEINIG, n. b. a. o. Prof. a. d. 
Techn. Hochschule Berlin, Aerodynamik der 
Luftschraube. XVI+ 484 S. m. 374 Textabb. 


ze. 1940, Verlag Julius Springer. Preis geb. 
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NACHRICHTEN 


Neu berechnete Tafeln. 


Auf S. 115 bis 129 seines Leitfadens zur Berech- 
nung von Schallvorgängen (Berlin 1939 bei Springer) 
gibt H. Stenzel Zahlentafeln für die Integrale von 
gewissen Funktionen, die mit den Besselschen 
Zylinderfunktionen Jp (x) für +p=2n+05 zu- 
sammenhängen (n=0,1,2,...). Es sei 

ine (— 1)" ae J+ (@n + 0,5) (x) 
08 


sin e 
= 
mit 
2n 1 
! 
und es sei 
x x 
mt) dt. 
a 


Die untere Grenze a ist so zu wählen, daß sich der 
Ausdruck für £ nur aus Gliedern mit dem Faktor 
sin z oder cos x zusammensetzen läßt. Die Tafeln 
geben die Werte der Funktionen o,, (2) und {sn (2) 
für n=0, 1, ....9, 10 und für z= 1,00; 1,25; 
...9,75; 10,00 mit 5 geltenden Stellen. 

Wenn mit allen 5 Dezimalen interpoliert werden 
soll, müssen noch die vierten Differenzen heran- 
gezogen werden. (Für lineare Interpolation dürfte 
der Schritt höchstens = 0,006 sein, für quadratische 
Interpolation höchstens = 0,04.) Beim Schritt 0,25 
kann der durch die Beschränkung auf die lineare 
Interpolation entstehende Fehler bis auf 0,008 an- 
wachsen, also schon die zweite Dezimale um 1 
fälschen. Bei quadratischer Interpolation kann der 
Fehler 0,001 erreichen, so daß man drei brauch- 
bare Dezimalen erhält. 
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Neu berechnete Tafeln. 


In seiner Züricher Dissertation „Über die Festig- 
keit der Kegelschale“ (Zürich 1917 bei Orell Füssli) 
setzt Fr. Dubois 


—il,(e), 
@y—ird)= 


und gibt auf S. 126/7 die Werte von /„, dI/dz, 
(/„dz mit 4 bis 9 geltenden Stellen für 2=0,6; 
1,1; 1,6; ... 38,5; 42,0; 45,5 (Schritt unregelmäßig). 
Auf 3 Seiten werden Kurven dieser Funktionen 
gegeben. 

[Bei F. Tölke, Besselsche und Hankelsche 
Zylinderfunktionen nullter bis dritter Ordnung 
vom Argument r Y i (Stuttgart 1936 bei Wittwer) 
findet man die Funktionen J, =— I,, Ja =—1,, 
vierstellig für r=2yx 
= 0,01; 0,02; ... 20,99; 21,00 (also z= 0,000 025 bis 
110,25). Aus diesen Zahlenfolgen kann man leicht 
die Ableitungen bilden und etwas mühsamer auch 
die Integrale.) 


Stuttgart. 
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Internationaler Mathematiker-Kongreß. 


Der internationale Mathematiker-Kongreß, der im 
September 1940 in Cambridge (Mass.) stattfinden 
sollte, ist auf einen günstigeren Zeitpunkt ver- 
schoben. 


Persönliches. 


Dem Direktor des Kaiser-Wilhelm-Institutes für 
Strömungsforschung in Göttingen, Prof. Dr. L. 
Prandtl, wurde anläßlich seines 65. Geburts- 
tages die Goethe-Medaille für Kunst und Wissen- 
schaft verliehen. 

Die Akademie der Wissenschaften in Wien hat 
das bisherige korrespondierende Mitglied, den o. 
Professor der Technischen Hochschule Graz, Dr. 
Karl Federhofer, zum ordentlichen Mitgliede 
ihrer mathematisch-naturwissenschaftlichen Klasse 
gewählt. 

Der o. Professor der Mechanik an der Technischen 
Hochschule Aachen, Dr. W. Müller, wurde zum 
o. Prof. an der Universität München ernannt, 

Es verstarb der frühere o. Professor der Geo- 
däsie an der Technischen Hochschule Berlin, Geh. 
Regierungsrat Dr. Karl Haußmann. 
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